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ABSTRACT 


An experimental investigation of a wing boundary layer subjected to periodic tur- 
bulent flow at a Reynolds number of 500,000 was conducted. Non-thrusting turbulence 
pulses were generated at a rate of 50 pulses per second with a turbulence intensity near 
10% bv spinning a rod in the free stream flow upwind of the test airfoil. 

Time-varying velocity measurements Were made at three representative chord lo- 
cations (laminar, transitional’turbulent, fully turbulent) at one angle of attack using 
single-element hot-wire anemometry. The characteristics of the boundary laver velocity 
profiles, turbulence intensity profiles, and velocity spectra and total spectral power were 
documented. 

The time-varving boundary layer response could be characterized by three flow re- 
gimes: undisturbed flow; turbulence pulse; and recovery period. The boundary layer 
exhibited a cyclic transition response varving from the undisturbed flow regime to the 
turbulence pulse regime back to the undisturbed flow regime. The turbulence pulse was 
found to penetrate into the entire boundary layer. The turbulence pulse prompted 
laminarization of the flow during the recovery period. Laminarization apparently re- 
sulted from rapid acceleration of the near-surface flow within the boundary laver due to 
turbulence induced momentum transfer down through the boundary layer and local flow 
acceleration following the Velocity deficit of the turbulence pulse. The effects of the 
periodic turbulence pulse were most noticeable in the recovery period at the 
transitional'turbulent and fullv turbulent regions of the wing. Recovery was generally 
characterized by lowered turbulence levels and reduction in boundarv layer thickness, 


indicating a stabilizing effect of the external turbulence pulse. 
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NOMENCLATURE 





c wing chord 
F tunnel calibration factor 
f. cutoff frequency 
f, Nyquist frequency 
f, samphng frequency 
K acceleration parameter —- ee 
tex 
p, total pressure 
p static pressure 
Ap differential static pressure 
q test section dynamic pressure 
R, chord Reynolds number —> 
, . u”? 
T, turbulence intensity — 
Us streamwise mean velocitv or ensemble average Velocity 
u-—s- streamwise instantaneous velocity 
v test section velocity 
X streamwise coordinate 
v _- coordinate direction normal to and away from surface 
6 boundary laver thickness 
v kinematic viscosity 
p density 
subscripts 
e boundary layer edge value 
co =o. free stream value 
superscripts 


fluctuating quantity 


time mean value 
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I. INTRODUCTION 


A. BACKGROUND 

Airfoil experimental research at low chord Reynolds numbers has become a major 
field of academic, military, and civil interest in the last few years due to the current de- 
velopment of low-speed remotely piloted vehicles (RPVs), manned and unmanned 
ultra-light aircraft, human-powered aircraft, and high altitude station keeping aircraft. 
Low Reynolds number flight (R, < 500,000) is generally considered as verv low speed 
flight, or flight with a short wing chord. Low Reynolds number flight also occurs at high 
altitudes where the air density is low. Most low Reynolds number data in the past have 
come from experiments on gliders and model airplanes. [Ref. 1] 

The aerodynamic response of a wing at low Revnolds numbers 1s considerably dif- 
ferent than for Reynolds numbers greater than about 500,000. At low Reynolds num- 
bers, laminar flow exists for a distance downstream from the leading edge. As such, the 
laminar to turbulent transition occurs much further back on the wing, often in the ad- 
verse pressure gradient region. The transition area is often quite broad and not well 
defined. In addition, the transition process is very sensitive to free stream disturbances, 
external turbulence, and surface imperfections. 

Laminar flow tends to separate easily at low Revnolds numbers especially in the 
presence of an adverse pressure gradient. Leading and trailing edge separation is fre- 
quently observed as well as separation along the airfoil. Separation along the airfoil 1s 
commonly referred to as the “laminar separation bubble”. Laminar separation bubbles 
greatly affect the performance of the wing. The most efficient low Reynolds number 
airfoils are designed to obtain maximum benefit from the separation bubble; but only 
because the separation bubble cannot be completely elinunated at low Reynolds num- 
bers. Numerous factors also affect the size, shape, and character of the separation 
bubble including free stream disturbances, vibration, surface condition, angle of attack, 
Etc. 

The critical nature of the boundary layer and the many factors-that influence it have 
prompted many recent investigations in the characteristics of airfoils and wings at low 
Reynolds numbers. Although investigations under steady state conditions provide im- 
portant lift, drag, and airfoil performance data, they may not reliably predict the laminar 


to turbulent transition, leading and trailing edge separation, and formation and reat- 


tachment of the separation bubble as these phenomena may be influenced to a high de- 
gree bv unsteadiness in the flow. [Ref. 1] 
Unsteady low Reynolds number aerodynamic research has recently been undertaken 


on three fronts: 
Airfoil response to high steady-state free stream turbulence 


Airfoil response to low to moderate frequency sinusoidal modulation of the free 
stream 


Airfoil response to transient free stream disturbances 


Meier and Kreplin [Ref. 2], Potter et al. [Ref. 3], and Kindelspire [Ref. 4], am@ng 
others, have investigated airfoil performance and boundary layer development and be- 
havior in the presence of high free stream turbulence. Meier and Kreplin found that 
both the intensity and spectral distribution (length scale) of the free stream turbulence 
affected airfoil performance. The most significant changes were noted when the length 
scales were on the order of the boundary layer thickness. Potter et al. also found that 
the free stream turbulence intensity spread into the entire boundary layer, but the effect 
was dependent on length scale. They also noted that the free stream turbulence has to 
penetrate and alter the velocity profile before producing anv significant change in the 
boundary laver detachment location. Kuindelspire found that increasing levels of free 
stream turbulence caused correspondingly earlier laminar to turbulent transition. He 
also confirmed that turbulence length scales on the order of the boundary layer thickness 
had the most influence on boundary layer thickness. 

Brendel [Ref. 5] performed a wind tunnel study with unsteady flow consisting of a 
streamwise sinusoidal velocity variation. Brendel found that, in general, the imposed 
unsteadiness had little effect on the mean boundary laver thickness. The unsteadiness 
did, however, reduce the thickness and length of the separation region associated with 
the separation bubble. He also found that lower frequencies (large length scales) had 
no effect on the overall flow. 

Howard [Ref. 6] most recently investigated the time varying characteristics of the 
boundary laver of an airfoil immersed in a propeller slipstream (see Figure 1). The 
propeller wake subjected the airfoil to a transient free stream disturbance. Howard 
found that the wing boundary layer varied from laminar to turbulent in response to the 
passage of the propeller wake. He also found, in general, that the effects of the propeller 
slipstream tended to enhance the stability of the boundarv layer and reduced the drag 


coefficient in the laminar portion of the slipstream cycle below its undisturbed value. 
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Figure 1. Wing Immersed in Propeller Slipstream: [Ref. 6: p. 17] 


In particular, Floward noted a distinct relaminarization process at work during the 
transition from turbulent to laminar flow as evidenced by a decrease in turbulence and 


an increase in velocity near the surface. 


B. THESIS OBJECTIVES 


The goal of this study is to thoroughly investigate the boundary layer response over 
a wing to periodic free stream disturbance at the laminar, transitional, and turbulent 
portions of the wing. 

This thesis describes the design, fabrication, and testing of a non-thrusting periodic 
turbulence generating device to simulate transient phenomena such as wind gust or tur- 
bulent wake disturbances. The investigation consists of wind tunnel tests of an airfoil 
operating in a periodically disturbed free stream at a chord Revnolds number of ap- 
proximately 500,000. Measured data are presented in the form of non-stationary mean 


velocity and turbulence time histories at various heights above the surface; as 


ensemble-average velocity profiles and turbulence profiles: and as velocity spectra and 
total power plots. 
These data will form a base for follow-on research at high angles of attack and for 


work within the separation region of low Reynolds number airfoils. 


II. EXPERIMENTAL EQUIPMENT 


A. APPARATUS 

The major hardware items utilized in this investigation included the wind tunnel, the 
Wing section and wing mount assembly, the periodic turbulence generator, the hot-wire 
anemometer, the traversing mechanism, and the data acquisition and processing system). 
This chapter presents the specifications, design, setup, and operating characteristics of 
all equipment used. 

1. Wind Tunnel 

The Naval Postgraduate School low speed horizontal wind tunnel was emploved 
for this experiment. Figure 2 shows the tunnel layout. 

The tunnel 1s of single return design and measures 64 feet in length, 25.5 feet in 
width, and supports a 32-inch by 45-inch by 4-foot test section. The tunnel was designed 
by Aerolab Development Company of Pasadena, California and has been in continuous 
use at the Navy School since the middle 1950's. [Ref. 7] 

The tunnel consists of five basic units: the fan section; diffusers; settling cham- 
ber; contraction cone; and test section. The tunnel 1s constructed of wood throughout 
and recently received a coat of fire resistant paint on the interior surface to meet OSHA 
standards. The paint roughened the surface somewhat; however, the level of free-stream 
turbulence was not significantly worsened due to the numerous anti-turbulence features 
built into the tunnel. 

The tunnel 1s powered by a 100-horsepower electric motor driving a three-bladed 
variable-pitch propeller. Tunnel velocity is controlled by electrically setting blade pitch 
angle at the control unit located adjacent to the test section. A four-speed Dodge 
transmission 1s installed between the motor and the fan to provide greater tunnel velocity 
control over the design range of the tunnel (up to 290 ft/sec). The transmission was left 
in 3rd gear for all data runs with a fan shaft speed of 750 RPM. 

Eight radial stator blades are mounted directly downstream of the fan on the 
drive shaft housing. These stator blades serve to straighten the flow and remove the 
swirl generated by the fan. These are the first of the installed anti-turbulence devices. 

The diffuser sections are located on either side of the fan section. The gradual 


increase in cross-sectional area and the long length of the diffuser sections gradually slow 
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Figure 2. NPS Low Speed Wind Tunnel 


the air while introducing little additional turbulence and Keeping friction losses to a 
muninium. 

The large settling chamber, located at the left of Figure 2, damps any remaining 
flow unsteadiness. In addition, two fine wire screens are mounted s1x inches apart in the 
settling chamber. These serve to further break up turbulent fluctuations in the flow 
prior to the air entering the test section. 

The contraction cone accelerates the air and delivers it to the test section. The 
large contraction ration of 10:1 helps assure low turbulence in the test section. 

The test section of the low speed wind tunnel is equipped with a full size re- 
flection plane mounted 4 inches above the tunnel floor. The reflection plane 1s perma- 
nently mounted and has an electrically-driven turntable mechanism 1n the center (not 
used in this investigation). The test section is additionally equipped with full length 
lighting located in the corner fillets and has a hinged window on either side. The walls 


of the test section are slightly divergent to counteract the contraction resulting from 


boundary layer growth. Breather slots are installed immediately downstream of the test 
section to maintain the test section at atmospheric pressure and to allow the tunnel to 
breath as necessary on start up and shutdown. 

Two sets of static pressure port rings are installed in the tunnel and used to 
measure the tunnel speed. One ring set is installed just upstream of the contraction cone 
and the other set is installed at the test section entrance. The pressure differential across 
the two ring sets is measured on a water micro-manometer. This pressure differential, 
once calibrated, corresponds to the velocity in the test section. The wind tunnel velocity 
calibration procedure 1s covered in Chapter III. 

2. Wing Section and Wing Mount Assembly 

The wing section and wing mount assembly are shown in Figure 3. This in- 
stallation was designed by Kindelspire [Ref. 4: pp. 11-15] and a complete description of 
the apparatus is contained in his work. The following paragraphs highlight the pertinent 
aspects of the wing and wing mount system. 

The airfoil-tested was a section of a Bell helicopter tail rotor (model and type 
unknown) having a 10-inch chord and 24-inch span. The airfoil shape 1s shown in Fig- 
ure 4 and the section coordinates are given in Appendix A. 

The airfoil section was stripped to bare metal, filled, sanded, and finished with 
high gloss black lacquer prior to installation. The lacquer finish was polished using fine 
polishing compound and given a final coat of automobile wax to bring it to the highest 
possible gloss. A glossy, reflective surface was important as the hot-wire anemometer 
probe was positioned to within 0.0032 inches of the surface by viewing the probe and its 
reflection through a sighting level. The probe positioning procedure is discussed in 
Chapter III. 

The wing mount assembly consists of two 8-in wide machined aluminum 
endplates mounted Vertically between the reflection plane and the top of the test section. 
The wing section 1s mounted between the endplates on rotatable inserts that allow the 
angle of attack (AOA) to be varied from +35 degrees to -15 degrees. Graduated 
markings on the rotatable inserts are used for coarse setting of AOA. The acutal AOA 
is determined by measuring the height of the leading and trailing edges above the re- 
flection plane and solving analvtically for AOA. 

All data collected during this experiment were taken l-inch left of mid-span 
(looking from the rear) due to the presence of several mars in the wing surface at the 
mid-span line which were caused by the measuring micrometer when the wing section 


coordinates were documented. The l-inch offset was sufficient to remain free of dis- 
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Figure 3. Wing and Wing Mount Assembly 


turbances generated by the surface imperfections! while also remaining clear of any end 
plate effects [Ref. 4: pp. 13-14]. The wing was set at 0 degrees angle-of-attack (AOA) 
on the inscribed end plate graduations. This resulted in an actual measured AOA of 0.71 
degrees. All data were collected at this AOA. 
3. Periodic Turbulence Generator 

Periodic turbulence pulses were generated by spinning a cylindrical metal rod in 
the freestream just upwind of the wing section. The turbulence generating assembly, 
which resembles a propeller, is shown in Figure 5. The mechanism consists of the base 
support assembly, cylindrical bar and hub, flexible drive shaft, and motor and 
timing’'position sensor. The support assembly and cylindrical bar’ hub were manufac- 


tured in the NPS Department of Aeronautics and Astronautics machine shop facilities. 


1 A general rule-of-thumb is that disturbance propagation 1s confined to a 15-degree arc 
downstream of a surface imperfection. 
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Figure 4. Experimental Airfoil Section 


The support assembly was constructed out of aluminum stock and designed to 
bolt to the reflection plane such that the hub was position 2.5 inches below and If inch 
ahead of the wing at center span. The hub was placed well below the wing in order to 
allow the disturbances generated by the hub section to pass under the wing. The leading 
edges of the mount were also rounded and/or beveled to minimize drag and free stream 
disturbances. Figure 6 shows the wing/wing mount and periodic turbulence generator 
installed in the NPS wind tunnel. 

The vertical strut of the support assembly was bored to accept the output end 
of the flexible drive shaft. The drive shaft was held firmly in place by a bracket mounted 
on the backside of the vertical strut. 

The bar and hub were both machined out of stainless steel. The 18-inch by 
3/8-inch diameter bar was secured to the hub bv a standard 1/8-inch roll pin. The hub 
was bored and tapped to screw onto the threaded end of the flexible drive shaft extend- 
ing from the support assembly. The hub was secured in place on the threaded shaft by 
a locknut. In addition, the counterclockwise direction of rotation was chosen such that 
the hub would tend to tighten during operation. 

The 38-inch flexible drive shaft used to couple the motor and the bar/hub was 
purchased at a local hardware store as an electric drill extension. As noted previously, 
the output end had a 3/8-inch threaded shaft onto which the hub assembly was mounted. 
The input end of the flexible shaft was made to fit a 1/4-inch drill chuck. Use of the 
flexible drive shaft, rather than direct drive, reduced the size, weight, and complexity of 
the support assembly as it did not have to support the heavy drive motor. 

The drive motor and position sensor are shown in Figure 7. The motor was 


mounted on the tunnel floor well clear of the wing environment at the downwind end 
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Figure 5. Periodic Turbulence Generating Assembly 


of the test section. The motor is a 1 4-horsepower permanent magnet variable speed 
motor. Motor speed could be controlled over the range 60-2500 RPM. An adapter with 
a 3,8-inch threaded output shaft was manufactured to fit the 1,2-inch motor shaft. A 
1'4-inch drill chuck was threaded on the adapter and used to drive the flexible shaft in 
this installation. A ferro-magnetic transducer was installed adjacent to the motor shaft 
to sense shaft speed and shaft position from a toothed steel disk mounted on the motor 
shaft. The toothed steel disk, seen in Figure 7, was modified by removing all but two 
of the original 30 teeth. Thus modified, an electrical pulse was generated twice each 
revolution corresponding to each passage of the bar in front of the wing section. This 
signal was used to drive an electronic RPM counter and provided timing information for 
the data acquisition system and data analysis software. 
4. Hot-Wire Anemometer 
Velocity and turbulence measurements were made using a Thermo System, Inc 


(TSI) IFA-100 Svstem Model 150 Constant Temperature Anemometer and Model 157 
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Figure 6. Wing/Wing Mount and Periodic Turbulence Generator Installed in NPS 


Wind Tunnel 


Signal Conditioner. Figure § shows the IFA-100 svstem installation. The IFA-100 


system is used with touch panel entrv of all ancmometer and signal conditioner settings 
and digital panel readouts. 

The Model 150 Anemometer module was operated in the standard bridge con- 
figuration with a standard 15-foot RG-58A U probe cable. This configuration provided 
Cmircimely low noise output signals with a tvpical frequency response of DC to 150 KI lz 
[Ref. S: pp. 8-3 to $-4]. These parameters were important in order to accurately measure 
the low level, high frequency turbulence components of the velocity signal. 

The hot-wire signal was processed using the Model 157 Signal Conditioner to 
improve the velocity resolution of the measurement system and to provide low pass fil- 
tering. The unconditioned hot-wire output voltage over the calibration velocity range 


varied from 1.5 volts at a flow rate of 40 f{t’sec to 1.8 volts at 140 ft sec. With the given 


analog to digital voltage resolution of 2.44 millivolts, the resulting velocity resolution 
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Fieure 7. Drive Motor and Position Sensor 


was, at best. 0.8 ft sec. Bv applving a -1.0 volt offset and gain of 5 to the raw signal with 
the signal conditioner, the output voltage over the calibration range was made to varv 
from 2.5 to 4.0 volts giving a velocity resolution of 0.16 {[t sec--the improvement in re- 
solution being directly related to the gain setting. 

The low pass filter of the Model 157°Signale@onditioner was"set"as "elec eaae 
possible to the Nvquist cutoff frequency. The Nvquist*cutoff frequeney, {, ,1s*aeiied 


by the formula: 


f= (1) 


5 Fa 


Where f, 1s the digital sampling rate in Hertz (Hz). To avoid serious errors in the"@am 
It is important that the signal have no significant components in the frequency range 
above f, prior to digital conversion {Ref. 9: pp. 68-69]. Setting the low pass filter to f, 


effectively removes all the unwanted high frequency components. 
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Figure 8. IFA-100 Anemometer System 


A Dantec 55p11 straight general purpose hot-wire probe with a 5-micron di- 
ameter platinum plated tungsten wire was used for all measurements. 

The oscilliscope and digital voltmeter, also mounted on the IFA-100 equipment 
rack (Figure 8), were used to monitor the hot-wire signal during each run. The digital 


voltmeter was set to read the mean DC, or steady state, component of the signal and 


served to monitor the steady state velocity. The oscilliscope displaved the AC, or tur- 
bulent, component of the signal and served to monitor turbulence in the flow. 
The calibration procedures for the anemometer are covered in Chapter III of 
this report. 
5. Three Dimensional Traverser 

Figure 9 shows the Velmex 8300 Control/Driver 3-D traversing system em- 
ployed to position the hot-wire probe through the boundary layer. The traverser con- 
figuration used by this investigator was identical to that utilized and described by 
Kindelspire (Ref. 4: pp. 21-24]. | 

The traverser was mounted above the test section, thus allowing the hot-wire 
probe to be positioned spanwise, chordwise, and vertically in relation to the wing sec- 
tion. Minimum incremental movement in any of the three directions was 0.000125 
inches. The traverser could be controlled either manually or by software command from 
a personal computer interface. An IBM PC/AT was used to control the traverser in this 
installation. The software routine used to drive the traverser is explained in a later sec- 
tion of this chapter. 

6. Analog-to-Digital Data System 

Data collection was accomplished using a PC/AT-type personal computer cou- 
pled to a Metrabyte Model DAS-16F analog-to-digital converter. 

The PC/AT-type personal computer used for data collection and analvsis was 
configured with a 40 megabyte hard disk, two high density disk drives, and operated at 
a throughput speed of 12 MHz. The computer was additionally equipped with | mega- 
byte of random access memory. The large storage capacity was necessarv to handle the 
volume of data collected in this investigation. Each sample point yielded close to } 
megabyte of raw data and over 150 points were sampled. 

The Metrabyte DAS-16F is a 12-bit verv high speed analog to digital interface 
designed for use with IBM and IBM compatible personal computers. The board fits 
directly in an expansion slot in the computer. The analog inputs were brought in 
through an external interface junction box. The board operates at a maximum rate of 
100,000 samples’sec (Hz) with up to 16 single ended or 8 differential analog input 
channels. The DAS-1I6F offers a wide range of operating parameters including four 
triggering options, three data transfer options, nine selectable input voltage ranges, and 
various other input/output modes. 

The DAS-16F was configured to read two differential channels (hot-wire voltage 


and timing pulse) over the input range of 0-5 volts. This resulted in an analog-to-digital 
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Figure 9. 3-D Traverser 


resolution of 2.44 millivolts. The board was operated at its maximum sample rate of 


100.000 Hz (50.000 Hz per channel) during this experiment. 


B. SOFTWARE 

The computer plaved an important role in all phases of data collection and analvsis 
in this investigation. The entire experiment was conducted using two PC AT-tvpe 
computers. The software programs created for each phase of the investigation are de- 
tailed in this section. The computer code for each of the programs developed ty this 
author 1s contained in Appendix C through Appendix F. 

I. Data Collection Software 
The programs in this section facilitated data collection and performed the fol- 

lowing control functions: 

3-D Traverser Control (TRAVERSE.BAS) 

Hot-Wire Cahbration (HWCAL.BAS) 


Analog-to-Digital Board control (LABTECH NOTEBOOK) 
Data Acquisition in the Boundary Layer (BLDATA.BAS) 


a. Program TRAVERSE.BAS 

The program TRAVERSE.BAS provides discrete or automated motor 
movement control of the VELMEX §300 traverser through interactive keyboard com- 
mands from an IBM PC/AT computer. The TRAVERSE.BAS program written by 
Kindelspire [Ref. 4: pp. 27-30] was modified by this author to remove manual data entrv 
routines and to improve the accuracy of the motor stepping routines. 

Program TRAVERSE gives the operator the option to input individual 
discrete step commands from the keyboard or to elect automated movement sequences. 
Discrete inputs were generally used to position the hot-wire in the boundary laver prior 
to each investigation. Discrete movements could be entered for anv or all of the three 
directions with the minimum movement being 0.000125 inches in any one direction. 
Once the hot-wire was positioned, the operator then selected computer controlled 
movement to advance the probe through the boundarv layer during the data acquisition 
process. 

The computer stepped the hot-wire through the boundary laver according 


to the positioning algorithm of Equation 2: 


x)= BL} -[sin($) J} pie (2) 


Z(j) = Position in the boundary layer 


BL = Estimated boundarv laver thickness 
n = Number of data points desired 
} = Incremental data point 


The sine function based algorithm of Equation 2 is designed to concentrate 
many smaller incremental steps at the bottom of the boundary layer and gradually in- 
crease the increment size toward the top of the layer. Special care was taken in writing 
the actual stepping routine to ensure motor movement commands were given in multi- 
ples of the motor step interval (0.000125 inches) to preclude cumulative roundoff errors 
throughout the boundary layer. 

b. Program HWCAL.BAS 


This program, written by the author, is used to calibrate the hot-wire 
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anemometer. The program was written in Microsoft2 QuickBasic and compiled to gen- 
erate a stand alone executable program. 

The hot-wire is calibrated against a reference pitot-static probe located in 
the test section free stream adjacent to the hot-wire. Pitot-static input 1s in centimeters 
(cm) of water as read on a micro-manometer. These data are manually entered on the 
kevboard by the operator for each calibration point while the corresponding hot-wire 
voltage from the IFA-100 Anemometer is digitally read by the Metrabyte DAS-16F 
analog-to-digital converter board. The operator also enters the ambient temperature 
and pressure at the conclusion of eachrun. One thousand hot-wire samples are read and 
averaged at each calibration point3. The hot-wire was generally calibrated at 1, 2, 3, 4, 
6, 8, and 10 cm of water providing accurate calibration data over the velocity range 
40-140 ft/sec. Program output consists of five columns of data suitable for plotting us- 
ing most current graphing software packages: 

Pitot-static pressure (cm of water) 
Hot-wire voltage (volts) 
Free-stream velocity (ft/sec) 
(Hot-wire voltage) ** 2 


Melocitw.** .5 


The hot-wire calibration procedures are covered in Chapter III of this report. 
c. LABTECH NOTEBOOK a 

The commercial software package LABTECH NOTEBOOK (Version 4.3, 
marketed by Laboratory Technologies Corporation? provided set up and operating 
control of the DAS-16F analog-to-digital board. Using the NOTEBOOK menus, the 
operator was able to enter all the operating parameters including sample rate, channels 
to sample, run duration, start mode, output data file name and output device, and data 
output format. NOTEBOOK allowed this investigator to take and store up to 130,000 
data samples during any one run while operating at the 100,000 Hz rate. Other software 
control programs appeared to be limited to 32,000 samples under the same operating 
conditions. This feature was paramount, as a minimum data capability of 92,000 dis- 


crete samples was required at each point. 


2 Microsoft is a registered trademark of Microsoft Corporation 
3 The suggested analog to digital sample rate is 1000 Hz with a run period of 1.01 seconds. 


4 LABTECH 1s a registered trademark of Laboratory Technologies Corporation. 
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d. Program BLDATA.BAS 
This driver program controls the analog-to-digital data acquisition process. 
Program BLDATA.BAS was written in QuickBasic and compiled to generate a stand 
alone executable program. The user inputs the number of data points desired in the 
boundary laver and a common output file name. BLDATA then calls NOTEBOOK 
Which in turn initiates data collection. Once the data are collected and the output file 
written bv NOTEBOOK, BLDATA renames the output file giving it a sequentially 
numbered tag and prepares for another data point in the boundary layer. At the com- 
pletion of the run, the investigator has a directorv of sequentially named files relating to 
each data point. 
2. Data Analysis Software 
The author used the data analysis software routines specifically developed by 
Johnson [Ref. 10] to analyze the periodic, non-steady velocity and turbulence data ob- 
tained in this investigation. These programs are not reproduced here due to their length 
and considerable documentation. A summary, however, of Johnson's main programs 
as well as those written by the author 1s provided. 
a. Program CONVERT.FOR 
This program was written by the author to convert the raw hot-wire voltage 
and position sensor voltage into velocity data and a s¥ynchronization pulse. Input data 
was in ASCII integer form with the integer values 0 to 4095 representing 0 to 5 volts. 
The bar sensor voltage was processed to return a 2 at the exact center of the position 
sensor pulse and a 0 elsewhere. The hot-wire voltage to velocity conversion constants 
are obtained from the hot-wire calibration data from program HWCAL discussed pre- 
viously. 
b. Program MEANSMU.FOR 
This program, written by Johnson [Ref. 10: pp. 158-166], computes the time 
varving mean of a non-stationary, periodic time series. The data are first ensemble av- 
eraged in the time domain and then transformed into the frequency domain via a Fast 
Fournier Transform (FFT). The ensemble is then smoothed by applying a low pass dig- 
ital filter to the ensemble. The ensemble is then transformed back into the time domain 
by an inverse FFT and written to the output file. The amount of smoothing 1s a func- 
tion of the low pass filter cutoff frequency. The lower the cutoff frequency, the greater 
the smoothing and vice versa. The amount of smoothing can be specified by the user 
through keyboard input. The methodology of selecting the low pass cutoff frequency 1s 
detailed in Chapter IV. 
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c. Program TURBIN.FOR 
Written by Johnson [Ref. 10: pp. 173-181], this program computes the in- 
stantaneous turbulence intensity of a non-stationary, periodic time series. It is the tur- 
bulence intensity counterpart of the non-stationary mean program MEANSMU. 
TURBIN first calculates the time varving mean velocity of each ensemble by low pass 
filtering in the frequency domain. TURBIN then determines the instantaneous turbulent 
(fluctuating) component at each point on each ensemble. The turbulence components 
are then ensemble averaged and and written to the output file. As in MEANSMU, the 
user can specify the cutoff frequency of the low pass filter used to smooth the individual 
ensembles. 
d. Program PSD.FOR 
This program, written by Johnson [Ref. 10: pp. 182-190], computes the 
power spectral density of a real discrete steady time series. This program was for spec- 
tral analysis of the free stream flow in the test section. 
e. Program ENSPSDAV.FOR 
This program, written by Johnson [Ref: 10: pp. 191-202], performs spectral 
analvsis on a short time segment of an ensemble of non-stationary, periodic, turbulent 
velocity data. Whereas program PSD performs spectral analysis on a steady state signal, 
ENSPSDAY investigates the time varying spectrum on selected portions of the ensemble 
itself. As such, the author was able to separately investigate the undisturbed flow seg- 
ment, the turbulence pulse, and the laminar flow/recovery period and characterize the 


power spectrum and power content of each. 
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Hil. EXPERIMENTAL PROCEDURE 


A. WIND TUNNEL TEST SECTION CALIBRATION 

The purpose of wind tunnel test section calibration was to determine the relationship 
between the actual free stream dynamic pressure (q) in the test section and the measured 
static pressure differential (Ap) across the contraction cone. This relationship 1s used to 
indirectly determine test section free stream velocitv by measuring Ap. 

For this calibration, a reference pitot-static probe was inserted in the center of the 
test section. A differential micro-manometer was connected across the total and static 
ports of the pitot-static probe to measure q. The relationship between the pitot-static 


probe differential pressure, q, and test section velocity 1s given by Equation 3: 


y’ 


p ; 
¢=(i-P) =o (3) 





p, = Total pressure measured by the pitot tube 
p = Static pressure in the test section 
q = Test section dynamic pressure 
p = Density of air 
v = Velocity in the test section 
The static pressure differential across the contraction cone, Ap, was measured by a 
second differential micro-manometer connected across the two static port ring assem- 
blies mounted in the wind tunnel walls just before and after the contraction cone. In 


most cases q 1s linearly related to Ap and can be written: 
q = F(Ap) (4) 


where F is a constant called the Tunnel Calibration Factor [Ref. 7: pp. Al-] to Al-5]. 
F is determined experimentally by plotting the measured values of q and Ap. The slope 


F 1S CiVensaienen spelmt by: 


a a = lunnele@orreetionsef aetor (>) 


Once F has been determined, test section velocity can be calculated from Ap by: 
Zap 
= (6) 
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Calibration measurements Were taken at numerous Velocities over the range of 40 - 
140 ft,sec. The two differential pressures were recorded at each calibration point. 
Figure 10 shows the plot of q versus Ap for this calibration. The relationship is ob- 
served to be linear and the tunnel correction factor is 1.126. This value differed signif- 
icantly from the calibration factor of 1.075 which has been in common use at NPS since 
1982. The author speculates that the original 1982 calibration was made prior to in- 
stallation of the reflection plane. Another possibility 1s that the new calibration factor 


reflects the numerous small changes made to the tunnel over the years. 


B. HOT-WIRE ANEMOMETER CALIBRATION 

The purpose of hot-wire anemometer calibration was to determine the equation re- 
lating hot-wire output voltage to velocity. This equation was used 1n the data conver- 
sion program CONVDAT.FOR described earlier in Chapter III]. The hot-wire was 


calibrated before and after each data run. Hot-wire calibration consists of three steps: 
IFA-100 system setup and initialization 
Calibration data acquisition 


Curve plotting and least squares fit 


Initialization procedures for the IFA-100 system followed the steps for routine op- 
eration outlined in the instruction manual (Ref. 8: pp. 3-19]. The setup parameters for 
the IFA-100 were as follows: 


Bridge Standard Bridge Configuration 
High Pass Filter Not used 

Low Pass Filter 20 KHz 

Voltage Offset -1 volt 

Voltage Gain 5 


As mentined in Chapter IJ, the low pass filter was set to the Nyquist cutoff frequency 
calculated by Equation 1. The Nvquist cutoff frequency for this investigation was 25 
KHz. The low pass filter could only be set in increments of 10 KHz, however, so the 
filter was set to 20 KHz--the closest lower frequency. Voltage offset and gain were used 
to improve the measurement resolution of the analog-to-digital converter by expanding 
the output voltage range of the hot-wire anemometer, as discussed in Chapter II. 

The hot-wire was calibrated against a reference pitot-static probe located near the 
top of the test section. The hot-wire was traversed to a position adjacent to the pitot- 


static probe for each calibration. 
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Figure 10. Wind Tunnel Test Section Calibration 
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The hot-wire was calibrated at or near seven selected velocities (1, 2, 3, 4, 6. 8. and 
10 centimeters Ap across the contraction cone) using the software routine HWCAL.BAS 
described in Chapter IJ. The calibration velocities are centered about the expected op- 
erating speed of 95 ft,sec and extended over the range 40 - 140 {t/sec. 

The final calibration curve was obtained by plotting hot-wire voltage squared versus 
the square root of velocity. A nearly straight line relationship results for this fourth or- 
der fit. Most graphing routines will plot the data and provide the equation of the linear 


least squares best fit line. This linear equation is of the form: 
y = Ax + B (7) 


where y=./velocity and x = (hot-wire voltage)*. Using the constants A and B from 


Equation 7, the hot-wire calibration equation becomes: 
velocity = [A(voltage)’ +B} (8) 


This relationship can readily be used in any data conversion program such as 
DATCONYV.FOR (Chapter II) to calculate velocity. 


C. HOT-WIRE POSITIONING 

Great attention was given to the task of positioning the hot-wire in the boundary 
laver adjacent to the wing. The objective was to accurately position the hot-wire as close 
as possible to the wing surface while not damaging either the wing or the hot-wire. 

Vaught and Fernholz provide an excellent summarv of the various methods in use 
today [Ref. 11]. The optical sighting method was used for this investigation. Kindelspire 
also used this technique in prior work at NPS [Ref. 4: pp. 12-13]. 

Several modifications to Kindelspire’s methods were made by this investigator to 
improve positioning accuracy and repeatability. The techniques emploved are described 
below. : 

The height of the hot-wire above the wing surface was measured bv observing the 
hot-wire and its reflection in the wing through a surveyor’s sighting level. The sighting 
level was positioned approximately six feet behind and to the left of the tunnel test sec- 
tion. Figure 11 shows the sighting level in position adjacent to the NPS low speed wind 
tunnel. 

As noted previously in Chapter I], the wing was finished in high gloss black lacquer 
to enhance the reflected image of the hot-wire probe. The distance between the hot-wire 


probe and its image was determined by using the horizontal cross-hair and vertical offset 
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Figure 11.  Surveyor’s Sighting Level 


Vernier on the sighting level. The vertical vernier was adjustable over a range of + 45 
Vernier units from the center position. One vernier unit was determined to DC7Oigae2 
inches of vertical probe displacement using the traverser, with 0.000125 inches stepping 
accuracy, as a calibration reference. The sighting level could be repeatedly setae 
to | 2 vernier unit resulting in a positioning accuracy of 0.0016 inches. 

The actual probe to wing distance is 1.2 the distance observed between the twoumi- 
ages. The hot-wire images were brought to within three vernier units separation for data 
collected at the 70%> chord location resulting in an initial hot-wire distance of 0.0048 + 
0.0008 inches. The probe was brought to two vernier units separation for all other wing 
positions (50% and 30°0 chord) which equated to an initial distance of 0.0032 + 0.0008 
inches. 

An external flood light was directed at the wing section and the hot-wire probe 
during the sighting process to enhance the image in the wing. The polished high gloss 


black finish on the wing provided a clear and sharp reflected nnage comparable in qual- 


Hot-Wire Probe s==p> 


Reflection o=_p»> 





Figure 12. Hot-Wire Probe and Reflection 


Pemuomiite actual probe as shown i Figure 12. No correction was made for the shght 
angle of the wing surface due to camber as this angle was less than five degrees at the 


points examined resulting in a positioning error of less than three percent. 


D. ACQUISITION OF EXPERIMENTAL DATA 
1. Test Environment Conditions and Setup Parameters 
All test were conducted at a chord Reynolds number of 500,000. Tunnel ve- 


locity for a given Reynolds number is given by: 





Kew 
aC (9) 
R, = Chord Reynolds Number 
u = Viscosity 
po = Density 
é 


Chord Length 


Solving Equation 9 using typical NPS ambient conditions (temperature = 68 degrees F, 
pressure = 29.95 inches Hg) where R, = 500,000 and c= 10 inches yields a tunnel oper- 
ating velocity of 95 ft/sec. All data runs were made at this reference free stream velocity. 

The setup parameters that remained constant throughout the investigation are 


summarized below: 


Wing Angle of Attack (AOA) 

Wing Chord 

Free Stream Reference Velocity 
Chord Reynolds Number 
Turbulence Generator Motor Speed 
Frequency of Turbulence Pulses 


Period Between Turbulence Pulses 


0.71 degrees 
10 inches 
EO lit, Sec 
500,000 
1,500 RPM 
50 Hz 


20 milliseconds 


Free Stream Distance bel ween taulses 22.8 inches 

Analog Channels Sampled 2 channels 
Hot-Wire Voltage Channel 1 
Timing Data ; Channel 2 


Analog-to-Digital Rate 50 KHz per Channel 


Sample Period 0.92 seconds 


Number of Pulse Periods Sampled 45 
Number of Samples per Pulse period 1000 


Period Between Each Sample 20 microseconds 


2. Equipment Configuration 

Figure 13 shows the equipment configuration used for this investigation. As 
noted previously, two personal computers were used during the data acquisition phase. 
One computer was dedicated to controlling the traverser mechanism while the other 
managed data collection and processing. Two computers were necessary as the operat- 
ing programs and data (92,000 samples) entirely filled the free memory of the data ac- 
quisition computer. 

Data collection was quite routine once the hot-wire anemometer was set up and 
calibrated, and the hot-wire probe initially positioned in the boundary layer. From that 
point, the operator had only to monitor the operating environment (tunnel velocity, 
tunnel temperature, turbulence generator motor speed, and hot-wire output), step the 


hot-wire through the boundary laver using one computer, and collect data with the other 
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Figure 13. Equipment Configuration 


computer. Approximately seven minutes were needed to collect and store the data at 
any one position in the boundary layer. While the actual data sampling period was less 
than one second, the computer took over 5 minutes to transfer the data from random 
access Memory to permanent disk storage. 

3. Flow Characterization 

Data were collected to characterize the flow in the following regions of interest: 
Test Section Free Stream 
Periodic Turbulence Pulses 


Boundary Layer Response to Periodic Turbulence Pulses 


7d 


a. Test Section Free Stream 
Data were collected in the free stream of the test section with all equipment 
removed to determine the intensity and spectral content of the flow. This step served 
to document the steadv state baseline characteristics of the free stream flow at 95 ft/sec 
and also allowed this investigator to compare the free stream turbulence data obtained 
in this investigation with that obtained previously by other investigators. 
b. Periodic Turbulence Pulses 
Data were collected with just the periodic turbulence generator installed and 
operating in the test section. The tunnel and the turbulence generator were set at their 
normal operating parameters of 95 ft/sec and 50 Hz respectively. The hot-wire was po- 
siuoned at the approximate vertical position of the upper wing surface (approximately 
half way between the hub and the rod tip). Measurements were taken four, six, and 
eight inches downstream from the rotating rod. These three positions coincided with the 
30%, 50%, and 70% chord locations on the wing. These tests served to document the 
time-varving velocity, turbulence, and spectral characteristics of the turbulence pulse 
generated by the rotating bar. 
c. Boundary Layer Response to Periodic Turbulence Pulses 
Measurements were made in the boundary layer at the 30%, 50%, and 70% 
chord locations on the upper surface of the wing with the periodic turbulence generator 
operating at 50 Hz. The wing AOA was measured at 0.71 degrees and the free stream 
velocity was set at 95 ft'sec. Data were taken at up to 50 points over the estimated 
boundary layer thickness at each chord location. The spacing of the sample points taken 
at each of the three chord locations is given in Appendix B. 
All data measurements were taken over a four day period. One additional 
dav of tunnel operation was needed following the data collection runs in order to inves- 


tigate the hot-wire drift problems addressed in Chapter IV. 
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[Ive RESULTS 


A. HOT-WIRE ANEMOMETER CALIBRATION 

Hot-wire calibration was conducted before and after each boundary laver investi- 
gation to confirm the stability of the hot-wire anemometry system and to ensure the 
accuracy of the data. Roane [Ref. 12: pp. 40-41] experienced some calibration curve drift 
in previous hot-wire work in the NPS low speed wind tunnel. As such, a minor amount 
of drift was expected. This investigator planned to average the initial and final cali- 
bration curves and use the resulting average to convert each day's data. The tabulated 
data from the calibrations runs on each day tended to confirm the expectation of only 
slight drift and, therefore, the data acquisition runs were completed as planned. Subse- 
quent detailed analysis of the calibration data for day three, however, showed that the 
initial and final curves differed by approximately 14%. The drift for day two was less 
but still significant at approximately 8%. The drift for day one could not be determined 
as the final calibration data were lost. The results of the calibration runs for day three, 
presented in Figure 14, show the divergence between the initial and final calibration 
data. The magnitude of the error was too great to resolve by averaging. In addition, 
the cause of the drift and its characteristics were unknown. 

This investigator planned and conducted an additional six-hour tunnel run on the 
following day in order to document and characterize the observed hot-wire drift. The 
tunnel was run at 95 ft/sec with just the pitot-static and hot-wire probes installed in the 
test section. Calibrations were conducted at the start and everv hour thereafter. The 
results of this run, plotted in Figure 15, show the shift of the calibration curve over time. 
The drift rate is greatest at the beginning and becomes very small after six hours. The 
effects of the drift and the induced error are greatest at the higher velocities due to the 
fourth order relationship that exists between hot-wire voltage and velocity. 

The drift rate was graphically analyzed at a reference velocity of 100 ft/sec and ex- 
panded bv extrapolation to cover an eight-hour period. Figure 16 shows the change in 
calibration of the hot-wire over eight hours. Notice that nearly 50% of the drift occurs 


in the first two hours of operation. 


5 The final calibration data for day one was not saved due to operator error. 
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Figure 16. 


oe 


The author applied the drift vs. time relationship of Figure 16 to each davw’s cali- 
bration data in order to generate a set of time varying calibration equations. The 
equations were staggered such that the maximum velocity difference between sequential 
calibration curves was one ft’sec or less. These corrections assumed that the drift rate 
determined during the test run was representative of the drift experienced on each day 
and that the run duration for each days was eight hours. Corrections were applied only 
to the data taken on day two (50% chord) and day three (30% chord). 

This investigator believes the hot-wire drift is a result of environmental contam- 
ination of the hot-wire probe. Hot-wire contamination appears to be a well known 
phenomenon among experienced researchers, but a subject about which little 1s printed 
In current hot-wire operating manuals or in fluid flow measurement texts. An older 
publication by TSI has this to sav about hot-wire stability: 

For a mean flow measurement to be useful it must maintain its calibration over 
long time periods. Hot wires are notoriously troublesome in this regard due to their 
small size making them susceptible to picking up small contaminants .... The 
larger hot film sensor has been a great improvement in this regard. The effects of 
surface contamination tend to be a function of sensor size and a ceramic coating on 


the sensor surface prevents corrosion and erosion. Hence, the hot film is far supe- 
rior for mean measurements. (Ref. 13: p. 91] 


Figure 17, reprinted from the above TSI publication, shows the typical change in 
calibration of a hot-wire sensor and a hot-film sensor vs. time when subjected to an ar- 
tificially dirty flow stream. The drift vs. time curve for this investigation (Figure 16) is 
remarkably similar to the TSI drift curve. With this knowledge, it was not difficult to 
pinpoint the source of contamination to the wind tunnel itself. 

For the month prior to this investigation, the tunnel was used for a series of flow 
visualization studies using an oil-based smoke as a visualization medium. This type of 
medium tends to produce a certain amount of contaminants which, due to the closed 
circuit design of the tunnel, were trapped in the tunnel environment. This oil residue is 
the likely contaminant and cause of the observed drift. 

Little can be done about the tunnel contamination, but one likely solution to the 
drift problem 1s the use of a hot-film sensor instead of a hot-wire probe. As noted pre- 
viously, the hot-film is virtually impervious to particulate contamination. A second 
solution would be to pre-run the tunnel prior to taking data for six to eight hours to 


allow the hot-wire to stabilize. 
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Figure 17. Change in Calibration of Hot-Wire vs. Hot-Film When Subjected to a 
Dirty Flow Stream: [Ref. 13: p. 92] 


B. ANALYSIS OF TIME VARYING DATA 

Ensemble averaging methods were employed to characterize the non-stationary, pe- 
riodic data acquired during this investigation. The periodic nature of the data made it 
possible to collect a series of repetitive time histories, or pulse passage cycles, each as- 
sociated with the passage of one turbulence pulse. This series of time-histories 1s called 
an ensemble. A typical ensemble consisted of a series of 45 pulse passage cycles. The 
number of passage cycles in one ensemble was essentially limited by available computer 


memory. For this investigation, 45 cycles was the upper limit. Howard [Ref. 6: pp. 
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114-115] showed that 20 or more pulses were needed in each ensemble to accurately 
portray the actual mean values, substantiating that 45 pulses were sufficient in this case. 
1. Ensemble Characteristics 

Each of the 45 pulse passage cycles making up the typical ensemble for anv one 
data point is slightly different due to the random nature of the turbulent components 
of the signal. Figure 18 and Figure 19 show two sequential pulse passage cycles (in- 
stantaneous velocity vs. sample number) from a tvpical ensemble (Test Case A) taken 
near the wing surface at 70° chord. These data were collected using a single hot-wire 
sensor. The velocity data were sampled, or digitized, at SO KHz with the turbulence 
generator operating at 50 Hz. This digitizing rate results in approximately 1000 data 
points per passage cycle. The period of each passage cycle being 20 milliseconds, each 
sample represents 20 microseconds in time. Time is not recorded explicitly but can be 
calculated from the digitizing rate and sample number. The data are plotted with 
straight line segments between each of the sample points. 

Each passage cycle is typically made up of three representative flow regimes as 


depicted in Figure 18: 
Undisturbed Flow 
Turbulence Pulse 


Recovery Period 


The specific characteristics of these regimes are very dependent on chord location and 
position in the boundary layer; but the following general characteristics hold true for 
each. The undisturbed flow segment is representative of the natural steady-state flow for 
the wing when operated in an undisturbed free stream. Howard [Ref. 6: p. 81] confirmed 
this relationship experimentally by comparing the drag of the undisturbed segment with 
that for natural flow and finding the two equal. The turbulence pulse is clearly identified 
as a highly turbulent velocity deficit. The recovery period is broadly defined as the pe- 
riod from the end of the turbulence pulse to the point where the flow returns to the un- 
disturbed state. 
2. Mean Velocity Determination | 

The author employed program MEANSMU.FOR described in Chapter I] to 
extract the slowlv varying mean velocity from the ensemble of instantaneous velocities. 
In principle. MEANSMLU filters out the high frequency turbulence components while 


leaving the low frequency variations of the signal mean intact. The amount of smooth- 
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ing 1s determined by the user selectable cutoff frequency of the digital low pass filter in 
program MEANSML. 

Program MEANSMU extracts the mean through a two step process. First, the 
45 time histories are ensemble averaged to generate a single “representative” time his- 
tory. Ensemble averaging consists of simply stacking the 45 time histories one on an- 
other and taking the average of the 45 values at each sample point in the time history. 
Figure 20 presents the ensemble average of the 45 pulse passage cycles for Test Case 
A. Note that the large turbulent fluctuations observed in Figure 18 and Figure 19 have 
been reduced but not entirely eliminated. The three flow regimes are also quite pro- 
nounced and easily identifiable. 

Figure 21 through Figure 24 show the application of digital low pass filtering 
to the ensemble average of Figure 20 with low pass filter cutoff frequencies (f,) of 1082 
Hz, 541 Hz, 360 Hz, and 267 Hz respectively. Note that the amount of smoothing is 
directly related to f, ; the lower f, the greater the smoothing. 360 Hz was choosen as the 
optimal cutoff frequency based on the findings of Landrum and Macha [Ref. 14] in their 
analysis of the spectra documenting the growth of turbulent disturbances through the 
transition process. They found that the free stream energy 1s contained below 200 Hz 
while the turbulence energy is generally contained above 550 Hz in the early transition 
stages. The cutoff frequency was choosen to lie between these two frequencies. 
Coincidentiallv, the f.= 360Hz curve is relatively smooth and flat throughout the undis- 
turbed flow regime and displays good peak response over both the turbulence pulse and 
recovery period regimes. 

3. Turbulence Level Determination 

Turbulence values of the time-varying, periodic signals were calculated using the 
program TURBIN.FOR described in Chapter I]. TURBIN accomplishes this by calcu- 
lating the instantaneous turbulence at each point in each of the 45 time histories and 
then ensemble averaging the instantaneous turbulence magnitudes. The result is a 
time-varying representation of the mean turbulence at each sample point in the pulse 
period. It is important to note that TURBIN also employs user selectable digital low 
pass filtering in its routines. The author used 360 Hz in this routine to-maintain conti- 
nuity among the various output data. 

Figure 25 shows the turbulence trace for Test Case A. Note that the represen- 
tative flow regimes (undisturbed flow, turbulence pulse, and recovery period) are as 


clearly identifiable here as in the mean velocity traces in the previous section. 
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Figure 21. 
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Figure 23. 
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4. Spectral Analysis 

Spectral analvsis of the time-varving, periodic signals was conducted using pro- 
gram ENSPSDAV.FOR described in Chapter I]. ENSPSDAYV performs spectral analv- 
sis on a uSer selectable 256 sample-wide window of the pulse passage cvcle time history. 
In theory, spectral analysis can only be accomplished on a discrete (digitized) time series 
of infinite length. In practice, the window width, or number of samples, can be reduced 
through the application of special techniques to the Fast Fourier Transform--the com- 
promise being loss of frequency resolution with reduction in window width. The 256 
sample-wide window used in this investigation was selected as the best balance of these 
factors. Each flow regime fits into the window and the resulting frequency resolution is 
190 Hz. The reader must keep in mind that the output data from ENSPSDAV repres- 
ents the average taken over the entire window. Thus, dynamic changes that occur in the 
flow regime over the window period will be lost in the resulting average. 

Figure 20 and Figure 27 show the output from program ENSPSDAV for the 
turbulence pulse regime of Test Case A. Figure 26 displays the velocity spectrum of the 
signal. The vertical scale is a dB measurement of the relative turbulence energy with 0 
dB representing the highest magnitude component of the signal. The horizontal scale is 
frequency in Hz plotted out to the Nvquist frequency. Each 10 dB reduction in the ve- 
locity spectrum represents a 10 fold reduction in power level. For example, the power 
level at 4.500 Hz (-20 dB) is 1/100 the level of the highest magnitude component. The 
mean has been removed from the signal so that onlv the spectrum above 0 Hz 1s ana- 
Ivzed. I[t is important to note that the spectrum display provides no information about 
the absolute power contained in the spectrum, but only displays power levels relative to 
the highest magnitude in the spectrum. Figure 27 displavs the total (cumulative) power 
distribution plotted with regard to frequency. The total power curve asymptotically 
approaches the value of the total power in the spectrum. The shape of this curve is in- 
dicative of the power distribution within the spectrum. The sharp knee observed in this 
plot indicates that the low frequency components of the signal contain most of the 


spectral power with very little power in the higher frequency components. 


C. FREE STREAM ANALYSIS | 
The results of the free stream analysis are presented for flow at 95 ft/sec with all 


equipment removed from the test section. Figure 28 is a plot of velocity response over 


6 The mean is removed for all the velocity spectrum and total power distribution plots in this 
report. 
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Figure 26. Velocity Spectrum for Test Case A, Turbulence Pulse 


1000 samples. The signal is flat with a high frequency turbulent component of small 
ainplitude superimposed on the mean. The free stream turbulence was calculated at 
0.182% by the data analysis methods of this investigation. This figure compares favor- 
ably with the free stream turbulence intensity of 0.23% observed by Roane [Ref. 12: p. 
41} at a velocity of 116.4 ft/sec. The slightly lower value of free stream turbulence ob- 
served here 1s not unexpected due to the slower free stream velocity. 

Figure 29 shows the velocity spectrum of the free stream. The curve shows a sharp 
drop in intensity with increasing frequency indicating that the major frequency compo- 


nents are at the lower frequencies. 
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Figure 27. Total Power Distribution for Test Case A, Turbulence Pulse 


The total power distribution of the free strcam is depicted in Figure 30. The rapid 
rise and sharp knee of this curve confirm the concentration of power mm the lower fre- 
quencies. The frequency components above 2 KHz provide little additional power as 


shown by the flatness of the curve above this point. 


D. TURBULENCE PULSE ANALYSIS 

The turbulence pulse was examined at three positions (4, 6”, and 8”) downstream 
of the spinning rod with the wing and wingmount assembly removed. The smoothed 
mean Velocity tine histories for the three positions are shown in Figure 31. The plots 
show a nearlv symmetric velocity deficit trace at all three positions. Pulse duration 1s 


roughly 250 samples or 5 milliseconds and the deficit maxima are constant at -20 ft,'sec. 
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Figure 29, Free Stream Velocity Spectrum 


The width of the deficit appears to increase slightly with increasing distance downstream. 
peace locity deficit trace observed for the spinning rod is very similar to that generated 
by a stationary cylinder in steady cross flow. This similarity shows that spinning the rod 
does not significantly alter the classical response of the cylinder. 

Figure 321s the plot of the smoothed mean turbulence traces of the turbulence pulse 
at the same positions. The pulse is intensely turbulent with a peak Value of 10.5%. The 
turbulence traces coincide with the velocity traces; however, the turbulence traces be- 
come somewhat non-svmmetric downstream. In addition, the turbulence appears to 
anticipate or lead the velocity deficit by up to 60 samples or 1.2 milliseconds. Most 


notably, the turbulence magnitude decreases quite rapidly with distance downstream. 
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Figure 30. Free Stream Total Power Distribution 


This decay 1s not unexpected and 1s generally attributed to dissipation of the high fre- 
quency turbulence components due to viscosity [Ref. 15: p. 48]. 

Figure 33 shows the velocity spectra of the turbulence pulse regime and the undis- 
turbed flow regime taken at 4° downstream from the spinning rod. The turbulence pulse 
trace decays more rapidly and to a lower level than the undisturbed flow trace. Addi- 
tionally, the undisturbed flow trace levels out at approximately -20 dB for frequencies 
above 5 KHz while the turbulence pulse trace bottoms out at -42 dB. These plots show 
the turbulence pulse to be a low frequency dominated event while the undisturbed flow 


regime has components at much higher frequencies. 
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Figure 32. 
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Figure 33. Turbulence Pulse Velocity Spectrum 


The total power curves for the two regimes are plotted in Figure 34. The intensity 
of the turbulence pulse is clearly shown here; the pulse power is 4 orders of magnitude 
greater than the undisturbed flow. The rapid rise and sharp knee of the pulse trace also 


confirms the intense, low frequency nature of the turbulence pulse. Very little power 1s 
added by the frequency components above | KHz. 
FE. BOUNDARY LAYER ANALYSIS 


The response of the boundary layer to periodic turbulence was investigated at three 
locations on the upper surface of the wing: 30% chord, 50% chord, and 70% chord. 


These three locations were chosen as the natural steady-state (undisturbed) flow was 
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laminar at 30% chord, transitional/turbulent at 50% chord, and fully turbulent at 70% 
chord. 
1. 30% Chord 

The ensemble-average velocity and turbulence profiles are shown in Figure 35. 
These plots consist of a series of trme-sequenced profiles, or “snapshots”, taken at several 
points in the time history. These profiles depict the dynamic characteristics of the 
boundary laver throughout one pulse passage cycle. 

The profiles show that the flow at 30% chord transitions from undisturbed 
laminar flow with very little turbulence (sample 100); to highly turbulent flow with a 
nearly constant level of turbulence at the turbulence pulse center (sample 275); and back 
again to laminar flow (sample 550). It is interesting to note how quickly and uniformly 
the boundary laver responds to the turbulence pulse. The turbulence profiles show a 
nearly uniform turbulence increase and corresponding decrease with the passing of the 
turbulence pulse. At the same time, the velocity profiles transition quickly from laminar 
to turbulent and back to laminar in step with the turbulence pulse. The velocity increase 
observed near the surface of profile 350 1s unusual. The velocity increase occurs just 
after the turbulence pulse and is still evident in the following profile (sample 400). The 
possible reasons for this behavior will be examined later in this section. 

The smoothed mean Velocity time histories and the smoothed mean turbulence 
time histories of the pulse passage cycle are shown in Figure 36 and Figure 37 respec- 
tively. These traces represent the time histories at increasing vertical distances (v) above 
the wing surface. The bottom trace shows the flow history taken closest to the wing 
surface (near-surface region) while the top trace depicts external flow outside of the 
boundary layer. 

The velocity and turbulence traces in the external flow region are tvpical of the 
response observed downstream of the periodic turbulence generator (Figure 31 and 
Figure 32); 1.e., a symmetric, bell-shaped pulse followed by an immediate return to the 
free stream or undisturbed state. 

The velocity time historv traces show a Velocity deficit of 11 ft/sec in the ex- 
ternal flow due to the turbulence pulse. This deficit increases as the surface is ap- 
proached until very close to the wing. Quite near the surface, the velocity deficit 
decreases significantly and is followed by a sharp velocity increase just after the turbu- 
lence pulse. This velocity increase remains for some time and is slow to dissipate. The 


turbulence traces, on the other hand, show that the turbulence intensity of the pulse is 
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Figure 35. 
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Figure 36. Mean Velocity Time Histories in the Boundary Layer, 30% 
Chord: (Fs=50 KHz) 
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Figure 37. Mean Turbulence Time Histories in the Boundary Layer, 30% 
chord: (Fs=50 KHz) 
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felt uniformly throughout the boundarv laver with onlv a slight reduction in magnitude 
and stretching of the pulse very near the surface. 

Figure 38 through Figure 42 show the velocity spectrum and total power dis- 
tribution plots for the three major flow regimes (undisturbed flow, turbulence pulse, and 
recovery period) at points in the boundary layer. Recall that the velocity spectrum plot 
display is relative to the highest component in the spectrum while the total power dis- 
tribution plot depicts the cumulative actual power with regard to frequency. Recall also 
that the mean has been removed from the spectral analysis. Hence, power changes in 
the mean due to mean Velocity increases or deficits are not reflected in the spectral 
curves. One must view the mean velocity and turbulence time histories along with the 
profiles and the spectral plots as a group in order to understand the boundary laver be- 
havior. 

Notice how little the turbulence pulse affects the spectral distribution very near 
the surface in Figure 38. The velocity spectra for the undisturbed flow and recovery 
period regimes are verv similar and the total power distribution curves for the two re- 
gimes are comparable in shape and magnitude. In addition, there is no unusual recovery 
period in this region. The turbulence pulse comes and goes without affecting the spectral 
characteristics of the flow very near the surface. Higher in the boundarv laver, 
Figure 40 and Figure 41, the high frequency components become more significant in 
the undisturbed flow and recovery period regimes while the total power contained in 
both regimes drops sharply. This indicates more spectral power (turbulence) 1s con- 
tained in the boundary laver very near the surface than in the intermediate regions or 
even in the external flow region. Note: The author is unable to explain the unusual 
shapes of the undisturbed flow velocity spectra in Figure 4] and Figure 42 or the ap- 
parent increase in undisturbed flow total power in Figure 42. 

Together, the velocity spectra and total power distribution plots show that the 
flow recovers very quickly from the turbulence pulse and the pulse does not significantly 
affect the characteristics of the recovery period. In other words, the undisturbed flow 
and recovery period velocity spectra and total power distribution curves are virtually 
identical at all positions in the boundary layer. 

The author believes an understanding of the mechanics operating in the 
boundary laver at this chord location is essential in understanding the more complex 
responses observed at 50% and 70% chord. The discussion that follows is the author's 


interpretation of the operating mechanics. 
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Figure 38. 
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Figure 40. 
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Figure 41. 
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Figure 42. 
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The data show a normally laminar flow being uniformly disturbed bv a strong, 
highly turbulent pulse from the external flow. The boundary layer responds to the 
passing pulse bv momentarily transitioning to extremely turbulent flow. Once the pulse 
has passed, the boundary layer rapidly returns to low turbulence laminar flow. Note 
however, that the velocitv time history graph (Figure 36) shows that the mean velocity 
quite near the surface increases in response to the turbulence pulse. Furthermore, this 
velocity increase remains in evidence for a considerable period after the velocity profile 
returns to laminar. The cause of this velocity increase 1s important. 

The turbulence time histories show there 1s considerable turbulence and mixing 
throughout the entire boundary laver in response to the pulse. This seems to indicate 
that momentum transfer down through the boundary laver 1s possiblv responsible for the 
velocity increase near the surface. The velocity time histories seem to support this 
supposition with regard to observed characteristics of the velocity deficit. The growth 
of the velocity deficit as the surface 1s approached 1s consistent with momentuni transfer 
down through the. laver. 

The boundary layer region very near the wing surface cannot transfer momen- 
tum due to the solid wing boundary. Thus, the near-surface region absorbs all the 
transferred momentum as a Velocity increase. This is evidenced by the sharp decrease 
of the velocity deficit and rapid velocity increase, or flow acceleration, just following the 
pulse. These two phenomena are clearly observed in the bottom two velocity time his- 
tories. Finally, the velocity increase remains in the near-surface region due to the posi- 
tive velocity gradient (f=) in the boundary layer and is slowly dissipated by the 
viscous effects of the laminar layer. 

In summary, the velocity increase in the near-surface region of the boundary 
laver appears to be the result of a downward transfer of momentum within the boundary 
laver induced by the high turbulence of the turbulence pulse. 

2. 50% Chord 

Figure 43 shows the ensemble average velocity and turbulence profiles at 50% 
chord. The boundary layer response at this position is quite unlike the simple laminar- 
turbulent-laminar response observed at 30% chord. To begin with, the undisturbed flow 
at sample 150 is transitional;turbulent with 5% to 6% turbulence near the surface. The 
flow becomes extremely turbulent at sample 350 with a high level of turbulence across 
the entire profile. During the recovery period, the flow becomes what appears to be fullv 
laminar by sample 650; the turbulence is extremely low over the entire region and the 


boundary layer thickness has decreased from around 0.08 inches at sample 150 to 0.04 


65 


(#) ny (%) ny (%) ny (%) ny (3%) MJ (%) ny (%) ny 
olgegrZzaoo0osgegrZzaodogsgegrezaodsegrazeoaedosgsdysr+r €OAdlIsSare2O OEY &v 2 dO 


















00'0. 
7 
OlLO~ 
020 
2 oc'O 
-JOQLUINA 
oss OSL os9 OSS Ose OS¢ . Ost ajduos 
, ee ee ee POllad oan Pee MO} 4 sawibay 
Auanoaay ssue;Nqiny paqinysipun MO} 4 
(996/))) n (996/) 0 (96/3) 1 (296/33) 1 (aes/y) 1 (296/33) (986/43) 1 
Oct Olt OG OLOC! Olt OG OLoOCt OL OG OLOCt Ol! O& OL OLt Ol OG OLOLCL Ol! OG OLOLt Olt OB’ OL 
. Pra 11 00°0 
: : A : 7 0Z'0 
: oc’o 





semble Average Profiles, 50% Chord 


. En 


43 


Figure 


66 


inches. Samples 750 and 850 show the flow returning to the original 
transitional/turbulent state as evidenced bv increasing levels of turbulence in the near 
surface region. 

The smoothed velocity time histories (Figure 44) clearly show the influence of 
the velocity deficit over the entire boundary laver. The corresponding smoothed turbu- 
lence time histories (Figure 45), however, do not show any turbulence penetration into 
the near-surface region. Apparently the high level of natural near-surface turbulence 
overshadows the effects of the turbulence pulse. The apparent answer to this dilemma 
lies in the velocity spectra and total power distribution plots of the undisturbed flow re- 
gime and the turbulence pulse regime near the surface. 

Shown in Figure 46 are the velocity spectra and total power distribution curves 
taken near the bottom of the boundary layer. This figure corresponds to the near- 
surface region of the mean velocity and mean turbulence traces. The undisturbed flow 
total power trace is quite broad banded with significant power contributions from the 
frequency components up to 4 KHz due to the high levels of turbulence generated near 
the surface. The turbulence pulse, on the other hand, is a low frequency dominated 
phenomenon having no significant spectral power components above | KHz. Thus, 
while the turbulence intensity of the undisturbed flow and turbulence pulse are almost 
equal in magnitude, the spectral makeup of each regime varies greatly. The distinct na- 
ture of the turbulence pulse makes it felt throughout the turbulent near-surface lavers. 
The effect of the turbulence pulse on the undisturbed flow is clearly seen in this figure. 
The undisturbed flow spectrum 1s altered by the turbulence pulse such that the recoverv 
period spectrum resembles the turbulent pulse spectrum. 

Figure 47 through Figure 5] depict the changes in the velocity spectra and total 
power distribution at increasing vertical distances 1n the boundary laver. Near the sur- 
face (Figure 46), the turbulence pulse clearly alters the recovery period spectrum by re- 
moving the high frequency components contained in the undisturbed flow regime. In 
the middle of the boundary layer (Figure 49), the recovery period spectrum begins to 
resemble that of the undisturbed flow showing the decreasing impact of the turbulence 
pulse on the flow. Finally, in the external-flow region (Figure 51), the undisturbed flow 
and recovery spectra match, indicating that the effects of the turbulence pulse are purely 
transient. Of interest are the very large spikes in the spectra of the recovery period plots 
for y<0.1250 inches. The author speculates that these are Tollmien-Schlichting waves 


generated during the laminar to turbulent transition. 
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Figure 47. Velocity Spectra and Total Power Distribution, 50% Chord, y = 0.0083” 
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The total power distribution plots show trends similar to the velocity spectra 
curves. Near the surface (Figure 46), the turbulent power 1s quite high for all three flow 
regimes. Although the turbulent pulse does not affect the recovery period power level, 
the frequency distribution of the power 1s clearly altered by the turbulence pulse. The 
undisturbed flow curve rises slowly to the asymptote while the recovery period curve 
rises sharply and levels off very quickly. Near the middle of the boundary layer 
(Figure 48 and Figure 49), the recoverv period turbulence power is much lower than 
that of the free stream. This shows that there is a definite recovery process operating in 
this region as the flow is energized back to the free stream condition. In the external- 
flow region (Figure 51), the undisturbed flow and recovery period cumulative power are 
nearly equal indicating that the effects of the turbulence pulse are purely transitory. 

The reverse-transition to laminar flow (variously called reversion, inverse or re- 
verse transition, relaminarization, and laminarization) observed in Figure 43 1s identical 
to the phenomenon observed by Howard [Ref. 6] for a wing immersed in propeller wash. 
It appears that the turbulence pulse triggers some mechanism in the boundarv layer to 
laminarize the inherently turbulent boundary laver in this investigation as the propeller 
wash did in Howard’s work. 

Narasimha and Sreenivasan, in their comprehensive survey of relaminarizing 
flows, identified four major reversion mechanisms: reversion by dissipation; reversion in 
stably stratified flows; reversion in highly accelerated flows; and relaminarization by ro- 
tation. (igi 16) 

Laminarization by sudden acceleration of the flow in the near surface region 
appears to be the mechanism at work in this instance. Analysis of the turbulence pulse 
regime of Figure 44 reveals instantaneous accelerations of 8,000 to 10,000 ft/sec? are felt 
by the flow near the surface. The flow appears to be accelerated both by the turbulence 
pulse itself during recovery from the velocity deficit and by momentum transfer down 
through the boundary layer (see previous section). 

One of the most widely used parameters that have proven to predict the onset 


and course of laminarization of a turbulent boundarv laver 1s: 


pees) 
Uz(x) 4x 


(10) 


Where K is the acceleration parameter, U,(x) is the edge free stream velocity, v the 


x . nt ee : 
——— the instantaneous velocity derivative. K, as used in 


kinematic viscosity, and qs 
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Equation 10, is purely a free stream parameter that carries no information at all about 
the boundary laver itself. Numerous investigators have proposed that values of K 
greater than 2 x 10-° to 3.5 x 10-° will give rise to laminanzation in low speed 
incompressible flows. [Ref. 16] 

In this investigation, the accelerated flow is observed in the boundary layer and 
in the external free stream as well. For this reason, the author chose to evaluate K lo- 
cally within the boundary laver as a means of predicting laminarization. Applied in this 


manner, K becomes: 
fo eee ee (11) 


where LU(x,y) 1s the local ensemble average velocity at height v in the boundary laver. 

Values of the acceleration parameter, K, range from 2.2 x 10-° at the near- 
surface region to 1.7 x 10-° in the upper regions of the boundary laver at 50% chord 
using local velocity values. The value calculated near the surface is clearlv suflicient to 
induce laminarization of the turbulent boundary layer at this position. It 1s interesting 
to note that K 1s much higher at 30% chord where it remains fairly constant at 4.4 x 
10-° over the entire boundary layer region; but K does not play a sigmificant role at 30% 
chord other than hastening the return to naturally laminar flow. 

3. 70% Chord 

The ensemble average velocity and turbulence profiles plotted in Figure 52 
show vet entirely different characteristics than observed at 30% and 50% chord. The 
undisturbed flow is fully turbulent with a high level of turbulence spanning the entire 
boundary laver region (sample 200). The turbulence pulse flattens the velocity profile 
while inducing high turbulence (6% - 8%) uniformly over the vertical region (samples 
350 and 410). During recovery, the velocity profiles become laminar like (sample 720) 
as the turbulence decreases in the middle and upper regions of the boundary laver. The 
turbulence traces appear transitional/turbulent in the recovery regime. Sample $50 
shows the flow becoming fully turbulent again with the boundary layer height increasing 


and with increasing turbulence in the lower and middle section of the boundary layer. 
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Note that the turbulence level near the surface remains high and fairly constant 
throughout the period. Clearly the laminarization process 1s at work here also, but the 
effects are greatly reduced. 

The smoothed velocity time histories plotted in Figure 53 show that the velocity 
deficit is felt throughout the boundary laver. The velocity deficit grows near the mid- 
layer and decreases again near the surface. Of interest in the velocity time histories 1s 
the nature of the velocity increase following the turbulence pulse. The velocity increase 
is much higher at 70% > chord than at the previous wing locations. Note also the long 
recovery period and slow rise and decay time associated with the velocity increase. 

Figure 54 shows the smoothed turbulence time histories. Notice the high tur- 
bulence levels that persist in the near-surface regions during the entire period. The trace 
at y = 0.0702 inches would tend to indicate that the turbulence pulse affects the mid- 
layer region more than the near-surface region at this chord location. 

The velocity spectra and total power distribution graphs are shown in 
Figure 55 through Figure 60. In the near-surface region (Figure 55 through 
Figure 57), the turbulence pulse mainly affects the velocity spectrum of the recovery 
period while total power remains high and nearlv equal over the three flow regimes. 
Toward the middle of the layer (Figure 58 and Figure 59), the total power levels of the 
undisturbed flow and recovery period regimes drop and the velocity spectra align more 
closely. It is interesting to note that the total power of the recovery period is less than 
that for the undisturbed flow in the middle boundary laver region. This indicates that 
there 1s a definte recovery process at work here. In the external-flow region 
(Figure 60), the turbulence pulse does not affect the velocity spectrum of the recoverv 
period. Here the power curves for the undisturbed flow and recovery period regimes are 
verv low magnitude and equal. 

The local acceleration parameter (K) from Equation 11, varies from 1.3 x 10°° 
near the surface to 1.9 x 10-6 at the middle of the boundary layer. The low value of K 
near the surface indicates that this region should not laminarize. This is borne out by 
the data; the near-surface region 1s dominated by high turbulence and the mean velocity 
responds very slowlv to the turbulence pulse. Conversely, the middle boundary laver 
region, With a value of K close to that necessary for laminarization, has a laminar-like 
response: a sharp drop in turbulence following the turbulence pulse with a corresponding 
sharp velocity increase. 

Apparently the turbulence pulse prompts the middle boundary laver toward 


laminar flow but the pulse is unable to fully penetrate into the thick, highly turbulent 
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Mean Velocity Time Histories in the Boundary Layer, 70% 
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near-surface region. Regardless, the turbulence pulse appears able to prompt a 


laminar-like response at this chord location, contrary to what one would intuitively ex- 


pect. 
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Vo eGONCLUSIONS 


The time-varying effects of periodic free stream turbulence have been investigated 
on a wing boundary laver at three representative chord locations (laminar, 
transitional turbulent, and turbulent) on the wing. An experimental method was devel- 
oped and tested to generate non-thrusting, uniform periodic turbulence pulses in the 
external flow over the test wing. Although the turbulence pulses were restricted to one 
test condition, 50 pulses per second and a turbulence intensitv of near 10%, they dem- 
onstrated the dynamic character of the boundary layer subjected to single-event turbu- 
lence. The results of this investigation are detailed below. 

The turbulence pulse was a short (5 millisecond), high intensitv (10.5% Tu), low 
frequency dominated (no significant spectral components above ] KHz), non-thrusting 
event. 

The turbulence pulse was shown to fully penetrate the boundary laver at all three 
chord locations. It was observed that the turbulence pulse suppressed the high fre- 
quency spectral components in the recovery period regime. Actual penetration was less 
near the surface in the turbulent flow region of the wing, but the effects of the pulse were 
Siecquite evident. 

The boundary layer time-response showed a distinct transition from undisturbed 
(laminar, transitional, or turbulent) flow to highly turbulent flow when subjected to the 
turbulence pulse at all three chord locations. 

The recovery process leading back to undisturbed flow was highly dependent on the 


nature of the undisturbed flow as shown below: 


Laminar: Quick boundary laver response with a rapid return to undisturbed laminar 
flow. 


Transitional’ Turbulent: Slower boundarv layer response than in the laminar region. 
Laminarization was evident during the recovery process followed by a slow return to 
undisturbed transitional'turbulent flow. The boundary laver thickness was decreased 
during the recovery period, and the turbulence intensity dropped to 1/6 of its normal 
undisturbed value near the surface. 


Turbulent: Slow boundary laver response. Laminarization is quite evident in the mid 
and upper portions of the boundary laver and evident to a lesser degree near the sur- 
face. There was a marked decrease in boundary layer thickness during the recovery 
period in this region before returning to the undisturbed turbulent value. 
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Laminarization appears to be mainly prompted by rapid acceleration of the flow 
internal to the boundary layer. The internal flow is accelerated both by the turbulence 
pulse itself during recovery from the velocity deficit and by momentum transfer down 
through the boundary layer. The velocity time-history traces show that the small ve- 
locity deficit associated with the turbulence pulse in external flow is amplified in the 
middle boundary laver region due to momentum transfer down through the layer. The 
near-surface flow experiences a rapid velocity increase as a result of the momentum 
transfer. 

The total spectral power content of the flow regimes (undisturbed, flow, turbulence 
pulse, and recovery period) was more dependent on vertical position in the boundary 
laver than on chord location. Near the surface, the total power of all three regimes was 
very high and nearly equal even though the spectral distribution of the power varied 
considerably. The middle boundary layer region, on the other hand, is characterized by 
a high power turbulence pulse, medium power level undisturbed flow, and a low power 
recovery period. The very low power level observed in the recovery period shows that 
the turbulence pulse drives the flow to a very low turbulent energy state. Lamuinarization 
occurs as there 1s insufficient spectral power to support turbulent flow at this point. The 
flow slowly gains power and transitions back to the undisturbed flow regime as the 
steady state characteristics once again prevail. At the edge of the boundary laver, the 
power level of the turbulence pulse remains quite high while the power levels of the un- 


disturbed flow and recovery regimes are low and nearlv equal. 
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VI. RECOMMENDATIONS 


This investigation demonstrated the stabilizing nature of periodic turbulence on the 
transitional and turbulent regions of the wing. The study was done at a single low angle 
of attack at three representative chord locations (laminar, transitional, and turbulent). 
Further experimental work should be done at additional chord locations and at higher 
angles of attack to further study the stabilizing effect of periodic turbulence on attached 
boundary layers and on flow near separation. 

More research also needs to be undertaken to better understand the mechanism at 
work in the laminarization process. Future studies should investigate the individual 
pulse passage cycles in greater detail to determine if ensemble-averaging techniques ad- 
equately represent the highly dynanuc laminarization process, or if statistical methods 
or a single pulse passage cycle analysis would be more instructive. 

Data collection was extremely time consuming with the present configuration. The 


following modifications should be implemented for future investigations: 


Store and process all data in binary integer format rather than in ASCII real format. 
This will reduce the mass data storage requirement bv a factor of nine and also greatlv 
reduce data transfer and processing time. 


Expand the random access memory of the data acquisition computer so that the data 
collection and traverser control software routines can be managed by one personal 
computer. This will allow for full automation of the test sequence. 


Use a hot-film probe instead of a hot-wire probe in the NPS low speed wind tunnel 
to avoid contamination induced drift problems. 


Decrease the digitizing rate. This investigation has shown that little useful spectral 
information is contained in the frequencies above 10 KHz. As such, the sampling rate 
could easily be reduced to 20 KHz with no loss of meaningful data. This step alone 
will reduce miass storage by 60% with a comparable reduction in data handling and 
processing time. 


Future research should also be directed toward investigating the effect of periodic 
or single-event turbulence on an airfoil near stall or in post-stall flight with regard to loss 
of attached flow and loss of control. Questions to be answered include: 

What is the effect of the observed stabilizing phenomena on flow near separation? 


Can the results be worked to an advantage at high angles of attack with regard to 
aircraft supermaneuverability or controllability? This is important as augmented lift 
from dynamic-stall maneuvering will have limited application unless the pilot can 
maintain control to transition quickly from one maneuvering state to another. 
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Directly related areas of interest include continued studies of the effects of wind 
gusts, ship airwake turbulence, mussile wake, and aircraft wake on aircraft stability and 
controllability. Future studies of canard wake/airfoil interaction at high lift seem perti- 


nent, especially in relation to the continuing supermaneuverability studies at NPS. 
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APPENDIX B. DATA SAMPLE SPACING 
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APPENDIX C. PROGRAM TRAVERSE 


REM seslevievlesevesiodesiesioviesiesic PROGRAM TRAVERSE. BAS Jerdetesevededeslededestedevevedede dese desevede dedledovdedesodeclevle deste 
REM 
REM THIS PROGRAM, WRITTEN BY D. KINDELSPIRE AND MODIFIED BY R. RENOUD 
REM (SEP 88) CONTROLS THE VELMEX 8300 3-D TRAVERSER. THE INITIAL 
REM PORTION OF THE PROGRAM IS DEDICATED TO OPERATOR-SELECTABLE OPTIONS 
REM (I.E. MOTOR VELOCITY, MOTOR ACCELERATION, ETC. ) AVAILABLE FOR 
REM EACH MOTOR. DEFAULT MOTOR VALUES FOR EACH OPTION ARE PROVIDED 
REM IN THE PROGRAMMING AND THE OPERATOR CAN SELECT THE DEFAULTS FOR 
REM ALL MOTORS WITH ONE KEYSTROKE. THEW PROGRAM GIVES THE OPERATOR 
REM A CHOICE OF MANUAL INPUT (MEANING ONE MOTOR MOVEMENT FOR EACH 
REM INPUT WITH INPUT IN INCHES) OR COMPUTER-CONTROLLED MOVEMENT 
REM THROUGH THE BOUNDARY LAYER. THE SINE FUNCTION-BASED ALGORITHIM 
REM USED TO COMPUTE THE MOTOR MOVEMENTS THROUGH THE BOUNDARY LAYER 
REM GIVES MOVEMENT COMMANDS IN INCREMENTS OF THE MINIMUM MOTOR 
REM MOVEMENT STEP SIZE (0.000125 INCHES). 
REM 
REM seslesle eves deveslesleslosleviedioslovedloseslosicleslese desterlontie seve Fese ese see to se eee se sevesevesesesese eee ceseseae see 7. ee 
REM 

REM OPEN THE COM PORT AND INITIALIZE THE MOTOR SETTINGS 

OPEN ''coml: 1200,n,8,1,rs,cs,ds,cd. FOR RANDOM AS #1 

REM SET MOTOR DEFAULT VALUES 

DATA 1000,1000,1000,2,2,2,0. O0@Wa5 ,0-000125 ,0. 000125 

READ V1, V2; V3,6R1, R2, R35, Clee 2c s 

DIM Z(100), Z1(100), HWV(100), RMS(100), VEL(100), NDVEL( 100), VRMS( 100) 
PI = 3.1415926544 

PRINT 

PRINT TT rectes fantentectonteaten!e eslectecte lon sleslecloclocloctectelocton' wtorle wleste cloctoclontonte slevecles 'eglole sfevles ale le oie ft 

PRINT ' * USER MUST SELECT "CAPS LOCK’ FUNCTION weal! 

PRINT oe acleslecleveslestectes favioslestoutestes! levleslestes ote eate ele vlestes ae elewotouton aafoutantaata tf 

REM DISPLAY MOTOR DEFAULT “SETTINGS 

PRINT i Wedevederlecleclevevedeslevevledetesedeseslotc ve deve declededevevlevededeslevede! 

PRINT " INITIALIZED VALUES FOR ALL MOTOR SETTINGS:" 

PRINT '"' VELOCITY = 1000 STEPS/SEC"' 

PRINT " RAMP(MOTOR ACCERBRATION) = 2 (6000 STEPS/SEC 
PRINT " rien INCREMENTAL UNTES ARE = 

PRINT " sestevesledeveslevicvlevedeslcclestesievlestesicviesicsiovesleslesteveslestecie clesleviesieslesc 

PRINT 

PRINT "DO YOU WANT TO INPUT YOUR OWN INCREMENT EACH TIME OR HAVE THE” 
PRINT "COMPUTER COMPUTE MOVEMENTS THROUGH A BOUNDARY LAYER FOR YOU 2?" 
PRINT 

PRINT tt Jedevestedevecteslewe ve NOTE werlesleslesie de glecte stevlesteslente tl 

PRINT FOR COMPUTED BOUNDARY LAYER MOVEMENT, YOU MUST KNOW: " 

PRINT ' 1. WHICH MOTOR WILL MOVE THE APPROPRIATE DIRECTION. ” 

PRINT " 2. HOW THICK (IN INCHES) IS THE B.L."' 

PRINT " 3. HOW MANY DATA COLLECTION POINTS DO YOU WANT IN THE B.L." 
PRINT " sedeskcaesevestcsedeslcvedcacsescslskesescslestcaeseseskeveslevese!! 

PRINT 

PRINT "NOTE!! USE MANUAL CONTROL TO INITIALIZE PROBE POSITION BEFORE" 
PRINT "' SELECTING COMPUTER CONTROLLED MOVEMENT. i 

PRINT 

INPUT "MANUAL CONTROL OR COMPUTER CONTROL (ENTER 'MAN' or 'CP')"; CONS 
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IF CONS = "CP" THEN 3130 
REM OPTION TO CHANGE DEFAULT SETTINGS OF VELOCITY OR ACCELERATION RAMP 
PRINT 


PRINT 

PRINT ' DO YOU WANT TO CHANGE THE VELOCITY OR ACCELERATION RAMP" 
PRINT " DEFAULT SETTINGS? (Y or N)" 

PRINT 

PRINT "IF 'NO', THIS PROGRAM WILL THEN LET YOU DEFINE THE" 


PRINT "DISTANCE YOU WANT TO MOVE (IN INCHES). IF ‘YES'," 


PRINT "YOU CAN CHANGE ANY OR ALL OF THE DEFAULT SETTINGS FOR ANY MOTOR. " 
PRINT 

PRINT 

PRINT 

INPUT "DO YOU WANT TO CHANGE ANY OF THE DEFAULT SETTINGS? (Y or N)"; DS 
IF DS = "Y" THEN 1520 

IF D$ = "N' THEN 2150 

REM 

REM = ee" OPERATOR SELECTED MOTOR VARIABLES *#%% 

PRINT 

PRINT 

INPUT "WHICH DEFAULT VALUE? (ENTER '1'FOR VELOC OR '2' FOR ACCEL RAMP)"; L 


ON L GOTO 1620, 1860 

PRINT "DO YOU WANT TO CHANGE THE DEFAULT VELOCITY? (Y OR N)" 

INPUT V$ 

IF V$ = "Y'' THEN 1620 

PRINT "DO YOU WANT TO CHANGE THE DEFAULT ACCELERATION RAMP? (Y or N)" 


IF R$ = "Y'' THEN 1920 

IF R$ = "N" THEN 1380 

PRINT 

PRINT 

INPUT "WHICH MOTOR VELOCITY DO YOU WISH TO CHANGE? (1,2, or 3)"; J 


ON J GOTO 1660, 1760, 1810 

PRINT 

PRINT 

INPUT "ENTER DESIRED VELOCITY OF MOTOR #1"; V1 
PRINT 

PRINT 

PRINT 

PRINT "DO YOU WANT TO CHANGE VELOCITY OF ANOTHER MOTOR? (Y OR N)'" 
INPUT VS 

DievS = Y THEN 1620 

IF VS = "N" THEN 1360 

PRINT ; 

PRINT 

INPUT "ENTER DESIRED VELOCITY OF MOTOR 2"; V2 
PRINT 

Gone 1710 

PRINT 

PRINT 

INPUT "ENTER DESIRED VELOCITY OF MOTOR #3"; V3 
PRINT 

GOTO 1710 

PRINT 

PRINT 

INPUT "WHICH MOTOR ACCEL RAMP DO YOU WANT TO CHANGE? (1, 2, or 3)"; K 
ON K GOTO 1900, 1990, 2050 


i 


PRINT 
PRINT 


INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #1"; 


PRINT 
PREND 


PRINT "DO YOU WANT TO CHANGE THE ACCEL RAMP OF ANOTHER MOTOR? (Y or N)?" 


INPUT RMS 
IF RMS 
IF RMS 
PRINT 
PRINT 


'Y" THEN 1860 
"N" THEN 1380 


INPUT "ENTER DESIRED ACCELERATION RAMP OF MOTOR #2"; 


PRINT 
PRINT 
GOTO 2350 
PRA 
PRINT 


INPUT "ENTER DESIRED. ACCELERATION RAMP OF MOTOR #3"; 


PRINT 

PRINT 

GOTO 1730 

REM 

REM DEFINE DISTANCE TO MOVE MOTOR 

PRINT 

PRINT 

PRINT 

REM INITIALIZE MOTOR INCREMENTS TO ZERO 


Series le wtewiantantsctontsntsnion! eve sleteslevevodedeves! eve ve dedevevesies! 
Werk DEFINE WHICH MOTOR YOU WANT TO MOVE 
vere NOTE! ! ! 
PRINT " **¥ 
PRINT " ** 
PRINT " ** 


weate 


weatent 


sdevedeslevesovesicsleveveveseslesesesvesesteslesicscsleslesicsiesicslestesesleslestvesesievesiese 


PRINT 


INPUT "WHICH MOTOR DO YOU WANT TO MOVE? (1,2, or 3)"; L 


eseseskes'es 


Rl 


R2 


R3 


A POSITIVE ('+') INCREMENT TO A MOTOR = *%*"" 
MOVES TRAVERSER AWAY FROM THAT PARTICULAR MOTOR 


-- MOTOR #1 MOVES THE PROBE UPSTREAM AGAINST THE FLOW 
-- MOTOR #2 MOVES THE PROBE TOWARD THE ACCESS WINDOW 9 **"' 
-- MOTOR 43 MOVES THE PROBE VERTICALLY DOWNWARD 


sosesveseseseseses 


cle vleveviesle deve leslie slo clonic slewlenlewlenecle clecleslentectewteslontectonlenieste 1 


wleste! 


oveeeevees cu 


ON L GOTO 2350, 2460, 2530 

PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #1?" 

PRINT '" avavdedededtdesese (ENTER DISTANCE IN INCHES) Seaesededesesesede |! 
INPUT I1 , 
PRINT 

PRINT 

PRINT "DO YOU WANT TO MOVE ANOTHER MOTOR ALSO? (Y or N)?" 
INPUT 

IF C$ = "Y" THEN 2300 

IF CS = "N" THEN 2600 
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2460 
2470 
2480 
2490 
2500 
Zon 
2520 
220 
2540 
2550 
2560 
2/0 
2580 
2) 0 
2600 
2610 
2620 
2630 
2640 
2650 
2660 
2670 
2680 
2690 
2700 
210 
Za 


PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #2?" 

REND | seessetesedecrse3t (ENTER DISTANCE IN INCHES) 2%? ededededede"! 

INPUT I2 

PRINT 

GOTO 2410 

PRINT 

PRINT 

PRINT "HOW FAR DO YOU WANT TO MOVE MOTOR #3?" 

PRINT WT Sevesededesesevede (ENTER DISTANCE IN INCHES ) seseresesedesedede" 

INPUT I3 

PRINT 

GOTO 2400 

PRINT 

PRINT 

REM DISPLAY OPERATOR SELECTED MOTOR VARIABLES 

PRINT !'! dededevesevedevedevedesedededevedededededesesedededededededese sede 

PRINT 

PRINT "SUMMARY OF OPERATOR INPUTS: " 

PRINT " MOTOR I VELOCITY = ": V1 

PRINT " ACCELERATION RAMP = "; Ri 

PRINT " INCREMENTAL DISTANCE = "; I1; ''INCHES" 
PRINT" . © MOTOR #2 VELOCITY ="; V2 

PRINT " ACCELERATION RAMP = "; R2 

PRINT " INCREMENTAL DISTANCE = "; 12; ‘'INCHES" 
PRINT " MOTOR #3 VELOCITY = "; V3 

PRINT " ACCELERATION RAMP = "; R3 

PRINT "' INCREMENTAL DISTANCE = "; 13; ''INCHES" 
PRINT 

PRINT '! vesedededevevedesesksededkaedesesededededesededededeseskedesededeve' 

PRINT 

PRINT 

PRINT "DO YOU WANT TO CHANGE ANY OF THESE VALUES? (Y or N)" 
PRINT 

PRINT "ENTER 'N' TO START MOTOR MOVEMENT. ENTER 'Y' TO RETURN" 


PRINT "TO VARIABLE SELECTION SUBROUTINE. " 

INPUT VS 

IF V$ = "Y" THEN 1360 

GOSUB 3060 

PRINT 

PRINT 

INPUT "DO YOU WANT TO INPUT ANOTHER MANUAL MOTOR MOVEMENT (Y or N)"; MS 
IF MS = "yy" THEN 2140 

PRINT 

PRINT DO YOU WANT TO INPUT COMPUTER CONTROLLED MOTOR MOVEMENT?" 

PRINT ' Sevevesedevesesesle NOTE! Pt dedvstvededestcdest " 

PRINT '" ALL PREVIOUS MOTOR INCREMENT INPUTS HAVE BEEN ZEROIZED." 

PRINT »PROGAM WILL LET YOU CHOOSE MANUAL OR CP-CONTROLLED B.L. MOVEMENT. " 
PRINT ‘sever? (IF 'NO', THE PROGRAM WILL END). ‘veove!! 

PRINT 

INPUT "DO YOU WANT COMPUTER CONTROLLED MOTOR MOVEMENT (Y or N)"; NS 

IF NS = "Y" THEN 1220 

PRINT 

PRINT 

PRINT 


oh) 


t? tout 


PRINT devevededesedededededededededededededededesededededesede|! 

PRINT " THE PROGRAM HAS ENDED. " 

PRINT ul devededevedevevevedeskevededesesesedevededesevesededevede |! 

END 

REM devvdedededs MOTOR MOVEMENT SUBROUTINE Jedededcdcdodedeste 

PRINT #1, "&': PRINT #1, EB: ‘Cl= "> Cie €2—= SeG2 ees ee 
PRINT #1, Ii=; 11; svVi—: a ERIE"; Rl; 

PRINT #1, 's12= 2°: V2="> V2 hoa 

PRINT #1, “I3="; 13; °:V3= 3; V3; :R3= 3 R3; 7G) 

RETURN 

REM Seaevesicataicscvedesedescacdede seseccdodeaerc se sesesescoces wae acer att aca sea aCe eee ee 

PRINT 

REM were COMPUTER CONTROLLED BOUNDARY LAYER MOVEMENT *%7eeectder 


PRINT 
PRINT 
PRINT " 
PRINT " 
PRINT " 
PRINT " weds 
PRINT " %%* DEFAULTED NEGATIVE MOTOR INCREMENT 
ERNE wk =6 (i.e. MOTOR #3 WILL MOVE UPWARD 
PRINT " BY ENTERING A (+) B.L. DISTANCE). 


PRINT " sevevcacdeskveveveserodededededededevevedesedevevededevevevesedevevevedesededede! 


PRINT 

REM SET INITIAL BL HEIGHT AND NUMBER OF DATA POINTS TO ZERO 
i = 0 

ili =O < 

PRINT 
INPUT 
PRINT 
INPUT 
PRINT 
INPUT 
PRINT 
INPUT 
ON MN 
PRINT 
PRINT 
REM MOTOR 1 CP-CONTROLLED MOTOR MOVEMENT 
ZC0)—s0 

i 

12 
13 
BOR J = slaroen 

Z(J) = BL *@@1 - SIN((CPI / 2) * (1 #19 - Diwe 
STP = CINTCG@ZQJ) - ZCI - 1)) 40012508. 0001S 

Z(J) = Z(J - 1) + STP 

I1 = -STP 
PRINT 
PRINT 
PRINT " 
PRINT 
PRINT " 
PRINT 
GOSUB 3060 
PRINT 


dedesedededcdededcdededededededcdcdededededededcdedededek descdescdedededededede|! 
deve NOTE !!! seve"! 
w%* COMPUTER CONTROLLED BOUNDARY LAYER 

MOVEMENT IS PROGRAMMED WITH A veel 


s: s'e t? 


seve!! 


seve! 


veut " 
5 weve 


"ENTER ESTIMATED INITIAL POINT DISTANCE ABOVE SURFACE ( INCHES)"; ZIP 
"WHAT IS THE B.L. THICKNESS (IN INCHES) THAT YOU WANT TO MEASURE"; BL 
"HOW MANY DATA POINTS DO YOU WANT TO SAMPLE IN THE B.L."; N 


"WHICH MOTOR WILL MOVE YOU THROUGH THE B.L. (1,2, or 3)"; MN 
GOTO 3380, 3730, 4090 


oe 
=o oy] 


Jerereseaedesededkedcdkededededesedesededededeicd | 


DATA POINT "; J 
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PRINT 


IF J = 


INPUT 
NEXT 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
IF RS 
IF RS 
PRINT 
PRINT 


N THEN 3620 

"PRESS 'ENTER' FOR NEXT MOVEMENT’; MOVES 
J 

"ALL MOVEMENTS COMPLETE" 


"YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER MOVEMENT (Y OR N)?" 


"IF 'Y', THE PROGRAM WILL TAKE YOU TO MANUAL CONTROL SUBROUTINE. "' 
"TF 'N', THE PROGRAM WILL END." 


"ANOTHER MOVEMENT"; RS 
"vy" THEN 1290 
"N" THEN 2990 


REM CP-CONTROLLED MOTOR MOVEMENT FOR MOTOR 2 
Z(0) = 0 


Jka 
1 
i 
BOR J 


Z(J) = 


STP = 


Z(J) = 


[2 = 

PRINT 
ERINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
GOSUB 
FeReN 1 
PRINT 


IF J = 


INPUT 
NEXT 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
IF RS 
IF R$ 
PRINT 
PRINT 
REM C 
Z(0) 

Ii = 


0 
0 
0 
=1TON 

eee le SINGOPT 2) (1 +03 - 1) & (CN - 1)))) 
GmemeezC tye) 2z(J = 1)) /..000125) * .000125 

Zee 1) 4 SIP 
-STP 


Setetevedededededcdededededededek sesesededetceye || 


tT 


DATA POINT "; J 


3060 


N THEN 3980 

"PRESS ENTER FOR NEXT MOVEMENT"; MOVES 
Jl 

"ALL MOVEMENTS COMPLETED" 


- 


"YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER MOVEMENT (Y OR N)?"; RS 


"IF 'Y', THE PROGRAM WILL TAKE YOU TO MANUAL CONTROL SUBROUTINE. "' 
"IF 'N', THE PROGRAM WILL END." 


"ANOTHER MOVEMENT"; RS 
‘Y THEN 1290 
"N'' THEN 2990 


P-CONTROLLED MOTOR MOVEMENT FOR MOTOR 3 
= 0 
0 


10] 


4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
4230 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 


13 
FOR J 
Zs) 

STP = 
ZT) 

I3 = 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
GOSUB 
PRINT 
PRINT 


IF J = 


ENBOT 


= 0 
= 0 


1 TON 
BL * (1 - SINC(PL./ 2) * Gi ts “a> (Ch eee 
CINT((Z(J) - ZCJ - 1)) / . 000125) =. COO 

Z(J - 1) + STP 
-STP 


uu seacsescacdedesdesdedesededesedesdescdedesedescse de! 


1 


DATA POINT "; J 


3060 


N THEN 4340 _ . 
"PRESS ‘ENTER’ FOR NEXT MOTOR MOVEMENT"; MOVES 


NEAT J 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 
IF RS 
IF RS 
END 


"ALL MOVEMENTS COMPLETED" 


"YOU WANT TO REPOSITION TRAVERSER FOR ANOTHER MOVEMENT (Y OR N)?"; RS 


"IF 'y', THE PROGRAM WILL TAKE YOU TO MANUAL CONTROL SUBROUTINE. " 
"IF 'N', THE PROGRAM WILL END." 


"ANOTHER MOVEMENT"; RS 
"y" THEN 1290 
"N" THEN 2990 
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REM Vevledestesleslevedevevesloveslesic sfesleste wtoute PROGRAM HWCAL. BAS Jeves sevleve stove le lewlewleslewlewloulecleulewleuls evevesles toute ulon ewe seve 
REM 
REM This program is used to calibrate the hot wire in the low speed wind 
REM tunnel. The hotwire is calibrated against a pitot/static probe 
REM located in the free stream adjacent to the hotwire. Pitot/static 
REM output (delta P) is in cm of water. Hotwire output from the IFA-100 
REM is digitally read by a Metrabyte 16f A to D board controlled by 
REM Labtech Notebook software. Suggested A to D speed is 1000 Hz with a 
REM period of 1.01 sec. The first 1000 data points are used by this 
Inuig program. | The program output file is directly compatable with 
REM Grapher” for plotting calibration curves and determining 
REM conversion equations. 
isle) 
REM tevedevevedestdevededevedevetede dete de nece de BeTe Fe ede TERE Te VOTE ETE TET TEETER TEE LETTER BERET B EES 
REM wevevevededevevevovede ceva toda dese decede Tete Fede Te PETE EE LONE TEE TE TERETE TE TE TE TERE ETE DE TE TESTED TE BETS TET TET TET CTE TOTS TOTS TTS 
REM 
REM VARIABLES: 
REM dp. : Pitot/static probe delta P (cm water) 
REM hwv : Hotwire voltage (volts) 
REM vfps: Tunnel velocity (feet per second) 
REM ofi: Output file name 
REM 
REM Wesvevevevededecovedodedevede Tere cede Terese te Feeds Tee Te TOTS Te ETE ETE ETE TE TST TS TS TEES TESTE TE TESTS TS OTS TS TS TOTS TEU TSS TS ESTOS 
DIM DP(20), HWV(20) 
R = 53.3: REM gas constant for air 
G = 32.174 
CLS : PRINT "Enter output file name’: INPUT ; OFIS 
OPEN OFIS FOR OUTPUT AS #2 
PRINT oe MT om water'", alin is "' yfas",,  Vipy . Be 
I= 1 
PRINT : PRINT "Enter cm of water for calibration point "; I 
Peer; DPC) 
PRINT : HWVCI) = 
REM Execute A to D conversion of hotwire using Labtech Notebook "go’ 
SHELL © nb go' 
REM Labtech Notebook writes output file to data. prn 
OPEN "c: nb data data.prn” FOR INPUT AS #1 
FOR J = 1 TO 1000 
INPUT #1, HWV 
HWV(I) = HWV(I) + HWV 
NEXT J 
CLOSE #1 
HWV(I) = HWV( I) / 1000 
CLS : PRINT © ‘hvw for point "; I; " ="; HWV(I) 
INPUT "another calibration point (y/n)"'; C$ 
IF CS = ''n" THEN GOTO 1250 
IT=I+1 
GOTO 1070 
CLS : INPUT "Enter ambient temperature (Deg F)"; T 
T = T + 460 


APPENDIX D. PROGRAM HWCAL 
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1270 
1 ie 10, 
1Z90 
1300 
1310 
L320 
1550 
1340 
150 
1360 
1370 
1380 
1350 


INPUT "Enter barometric pressure (in. Hg)"; P 
PRINT "Wait - Writing Output Data File’ 

P= P / 29.92 4) eee 

RO =P / (R* T * G) 

FOR1J. °= gia lO ust 

VFPS = SQR(2 * DP(J) * 2.048 / RO) 

PRINT #2, DP(J), HWV(J), VFPS, HWV(J) 2, SQR(VFPS) 
NEXT J 

CLOSE.#2 


CLS : INPUT "Another calibration run (y/n)"; C$ 
IF CS = "y'' THEN GOTO 1000 

PRINT "BYE" 

END 
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APPENDIX E. PROGRAM BLDATA 


REM 3 wvscedtdodevedevededeseweiss PROGRAM BLDATA. BAS Wvevededededededtdededededevede ade ce cede cede ce de ee Te eee eK 
REM 
REM This driver program controls the analog-to-digital data acquisition 
REM process using the software package LABTECH NOTEBOOK in conjunction 
REM with the Metrabyte 16f A to D card. 
REM This progarm allowes the user to enter the number of data points 
REM for the run. It initiates the data collection by running LABTECH 
REM NOTEBOOK program "go" directly from the nb directory via dos 
REM command. This control program also generates incremental names for 
REM the data files. 
REM 
REM  wevevevedevevededovededodedededk devevetedece ce Neve Lede FORE TERETE TE TT LC TCL TEPER, TOR TELE DETE DEDEDE LEPE LE TORE LCE TELE LOTTE 
CLS 

INPUT "Enter number of data points for this wing position (1-99)"; N 
pee = 1) OR {(N > 99) GOTO 1010 

REM *“*"* generating string segments for data file names 

BS = MIDS(STRS(1), 2): REM string number "1" 

ES = MIDS(STRS(N), 2): REM ending string number ‘n" 

pom) XXXXXX” 

EXS = ". DAT" 

CLS 

PRINT "Data files will be incremented from: "’ 

PRINT (X$ + BS + EXS); " To "; (X$ + ES + EXS) 

PRINT 

PRINT Enter common data file id (XXXXXX) for this run 

PRINT "(Six characters maximum - no extension’ 

INPUT F$ 

BELEN FS) = 6 OR LEN(FS)@< 1 GOTO 1120 

CLS 

PRINT N; "Data files will be generated and incremented as follows: " 
PRINT (FS + BS + EXS); " To "'; (FS + ES + EXS) 

PRINT 

INPUT ‘Are the number of points and file names ok (y n)"; CS 

IF CS = "n" GOTO 1000 

IF CS = "y' GOTO 1230 

GOTO 1190 

CLS 

FOR I = 1-TON 

Pie ie mPRINE ready for point "; I; © of "; N; ' points": PRINT 
PRINT "Insert New Data Disk in Drive A” 

PRINT "Strike any key when ready" 

BEEP 

CS = INKEYS: IF C$ = '" GOTO 1290 

REM **% generating incremented data file name 

IF (I > 10) OR (I = 10) THEN IS = MIDS(CSTRS(I), 2) 

eee OP THe N eis = (MEDSCSIKS(O), 2) + MIDSCSTRSCI), 2)) 

FI$ = (FS Loe © ) 


Pra 
SHE 


SHELL ' 


R1l$ = 


NT ' 
LL 


‘Doing A to D conversion and writing data file "; FIS$ 
cd. nb’: REM DOS command to ensure system in nb subdirectory 
pee) ': REM executing labtech notebook via DOS from NB subdir 
"ren A: DATA. PRN ' 
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1380 
oe, 0) 
1400 
1410 
1420 
1430 
1440 


R2$ = R1S + FIS 

SHELL R2$ 

NEXT I 

INPUT "Another run’; C$ 
IF cs = "y" GOTO 1000 
PRINT "BYE" 

END 
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APPENDIX F. PROGRAM CONVERT 
PROGRAM CONVERT 


Cede dededesede deve tedeve deve ve Fe de Fe Fede Fo HN Te Tee Tee TE Te Te ETE TE Te TE TET TENE TE TE TEE TS TET TE TE TE NE TE TE TENE TE Te FETE TC TET TE TEE TENCE 


Celgce ©) C2 C2 Conganca Ga C2 C2 


10 





THIS PROGRAM CONVERTS THE RAW HOT-WIRE VOLTAGE AND POSITION SENSOR 
VOLTAGE INTO VELOCITY DATA AN A SYNCHRONIZATION PULSE. INPUT DATA 
IS IN ASCII INTEGER FORM WITH THE INTEGER VALUES 0 TO 4095 REP- 
RESENTING 0 TO 5 VOLTS. THE BAR SENSOR VOLTAGE WAS PROCESSED TO 
RETURN A VALUE OF 2 AT THE EXACT CENTER OF THE POSTION SENSOR 
PULSE AND A O ELSEWHERE. THE HOT-WIRE VOLTAGE TO VELOCITY CONVER- 
SION CONSTANTS ARE OBTAINED FROM THE HOT-WIRE CALIBRATION DATA 
FROM PROGRAM HWCAL. BAS. 


a aute ant e's e's a vd r rs PY —s - a ae’ ae’ eae @) 2&2 ee 4 @, 20) 4 @'S &) 2 &D a = Pm a a a a Ge ae' a a 2e.'ae i a ae a's ae‘ aaa “a a =a 
eS OSS a Cae See eee eee eee seseseseseseseseseiesesoiesedesesedesesedesesesesesededsesvesesesesesetesodedestevetesetesesesteststsioleldeds 


-~ 
¢ 


REAL HWV, PROP, A, B, PPROP, TIC 
CHARACTER*12 INFILE, OUTFILE, CHOICE*2 
INTEGER 1,J,K 

DIMENSION DAT (50000,2) 

A=0 > 

B=0 


PRINT '(/A )',' ENTER THE INPUT DATA FILE NAME: ' 
READ (*,'(A)') INFILE 


PRINT '(/A )', ' ENTER THE OUTPUT DATA FILE NAME: ' 
ReaD) Cx,’ (A) SOUTFILE 


OPEN (10, FILE=INFILE, STATUS = 'OLD') 

OPEN (11, FILE=OUTFILE, STATUS = 'NEW') 

PRINT '(/A )', ' VEL = (A * HWV**2 + B) ** 2! 
1 


A acne" 
B t 


ees : B 
PRINT '(/A/)', ' CHANGE A OR B (Y/N) ? ' 
READ (*,'(A)') CHOICE 
ieGcnOnebe.EQ, N’)..GOTO.7 
Pea (7A), CUA 
READ (*,*) A 
ieee fA) , ' B= ' 
READ (*,*) B 


PPROP=0 
PRINT '(/A/)', ' READING DATA INTO ARRAY - WAIT' 
DO 10nd = 1, 50000 

READ (10, *, END=50, ERR=50) HWV, PROP 


DAT (J,1) = HWV 
DAT (J,2) = PROP 
CONTINUE 
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80 
30 


CONTINUE 
PRINT '(/A/)', ' DOING DATA CONVERSION - WAIT SOME MORE' 

Ho 20 T = eee 

DAT (1,1) = (A * (DAT (1,1) * 5. / 4095.) ** 2 + B) ** 2 
PROP = DAT (1,2) * 5. / 4095. 

TC = "0: 

IF (PROP .GT. 1.5 .AND. PROP .LT. 2. . AND. PROP .LT. PPROP) 


+TIC = 2. 


PPROP = PROP 
DAT (1,2) = TIC 
CONTINUE 


PRINT '(/A/)', ' WRITING OUTPUT FILE - ALMOST DONE! ' 
DO 30 1 = aS 

WRITE (11,80) DAT (i131), vAleG 

FORMAT (F11.4, F5.0) 

CONTINUE 


WRITE (*,*) ' # OF CONVERSIONS 9s : sige 
CLOSE (10) 

CLOSE (11) 

WRITE (*,*) ' ANOTHER RUN (Y/N) ?' 

READ (*,'(A)') CHOICE 

IF (CHOICE .EQ. ‘Y') GOTO 5 


END 
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